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ABSTRACT 


A demonstrated need exists for better wind field information over the open ocean, 
especially as a forcing function for ocean circulation models. Microwave scatterometry, 
as a means of remotely sensing surface wind information, developed in response to this 
requirement for a surface wind field with global coverage and improved spatial and 
temporal resolution. This development led to the 1978 deployment of the SEASAT 
Satellite Scatterometer (SASS). Evaluations of the three months of SEASAT data 
have established the consistency of SASS winds with high quality surface wind data 
from field experiments over limited areas and time periods. The directional ambiguity 
of the original SASS vectors has been removed by Atlas et al. (1987) for the entire 
data set, and the resulting SASS winds provide a unique set of scatterometer wind 
information for a global comparison with winds from conventional sources. 


A one-month (12 August to 9 September 1978) subset of these dealiased winds, 
in the western North Atlantic, is compared here with a conventional, pressure-derived 
wind field from the 6-hourly surface wind analyses of the Fleet Numerical Oceano- 
graphic Center (FNOC), Monterey, CA. Through an objective mapping procedure, the 
irregularly spaced SASS winds are regridded to a latitude-longitude grid, facilitating 
statistical comparisons with the regularly spaced FNOC wind vectors and wind stress 
curl calculations. The study includes qualitative comparisons to synoptic weather 
maps, calculations of field statistics and boxed mean differences; scatter plots of wind 
speed, direction, and standard deviation; statistical descriptions of the SASS-FNOC 
difference field, and wind stress curl calculations. 


The SASS and FNOC fields are consistent with each other in a broad statistical 
sense, with wide scatter of individual values about a pattern of general agreement. The 
FNOC wind variances are slightly smaller than the SASS values, reflecting smoothing 
on larger spatial scales than the SASS winds, and the SASS mean values tend to 
be slightly higher than the FNOC means, though the increase is frequently lost in 
the large scatter. Exceptions to the pattern of relatively small consistent variations 
between the two fields are the pronounced differences associated with extremely strong 
winds, especially during Hurricane Ella, which traveled up the East Coast of the 
United States during the latter part of the study period. These large differences are 
attributed mainly to differences in the inferred positions of the pressure centers and in 
the response at the highest wind speeds (> 20m/s). The large statistical differences 
between the SASS and FNOC fields, present under high wind conditions, may yield 
significantly different ocean forcing, especially when the strong winds persist over 
longer periods of time. Under less intense wind conditions, usually prevailing over the 
ocean, the two fields correspond well statistically and the ocean responses forced by 
each should be similar. 


Thesis supervisor: Dr. Carl Wunsch 
Cecil and Ida Green Professor of Physical Oceanography 
Secretary of the Navy Research Professor 
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1. Introduction 


The large-scale oceanic circulation has been extensively studied, especially over 
the past 100 years, and considerable advances have been made in theoretical under- 
standing of the processes involved. A critical part of this study is the attempt to 
understand the forcing that drives the circulation. The two major driving forces are 
the stress due to wind blowing over the ocean surface and the buoyancy-driven mo- 
tions due to thermohaline processes. The directly wind-driven circulation is confined 
to, and dominates, flow in the upper layers of the ocean. It is primarily horizon- 
tal motion, but convergence and divergence of this flow give rise to vertical motion, 
creating regions of upwelling or downwelling. Theoretical study of the wind-driven 
circulation was motivated initially by attempts to account for observed surface cur- 
rents, starting with the work of Nansen (1898), Ekman (1905), and Sverdrup (1947). 
Further progress was associated with attempts to explain westward intensification, 
reflected in strong western boundary currents (Stommel, 1948; Munk, 1950). This 
work established our basic understanding of how wind stress drives the ocean circu- 
lation, and subsequent developments are reviewed by Veronis (1981) and Pond and 
Pickard (1983). 

While wind-forcing is accepted as the primary driving force in the upper layers, 
themohaline forces also play a significant role in ocean circulation, especially in the 
lower layers where they tend to dominate. In this case, water movements are associ- 
ated with changes in density, either from temperature or salinity changes. In the ocean 
these changes normally occur as an increase in density at the surface, either directly 
through cooling or indirectly through freezing and subsequent production of more 
saline water. Such processes cause a buoyancy-driven vertical flow that is followed 
by horizontal flow due to continuity requirements or to geostrophic currents resulting 
from the changes in density. Thermohaline processes are frequently associated with 
changes in weather and climate, and may also be strongly influenced by the wind 
field, especially in terms of surface heat fluxes and downward mixing in the surface 


layer. Therefore, the wind field plays a most critical role, both directly and indirectly, 





in forcing the ocean circulation, and its specification is crucial to reaching a more 
complete understanding of that circulation through theory and model developments. 

However, the development and verification of these theories and their applica- 
tion to global circulation models are hampered by the current lack of good wind 
fields over the open ocean. The most frequently used fields for wind-forcing are av- 
eraged estimates of the wind stress, available both globally (Hellerman, 1967,1968; 
Han and Lee, 1983; Hellerman and Rosenstein, 1983) and for particular regions (for 
example: Bunker, 1976; Hastenrath and Lamb, 1977; O’Brien and Goldenberg, 1982). 
These climatological values are averaged across monthly, seasonal, and annual peri- 
ods, generally from all available ship reports in the chosen period, extending over a 
considerable length of time and employing some wind speed to stress conversion. In 
the case of the Hellerman and Rosenstein (1983) analyses, for example, ship reports 
from 1870 to 1976 were used to form monthly averages, and stress was computed 
from wind speed based on a wind speed and atmospheric stability dependent drag 
coefficient. Many uncertainties still exist in the details of the wind speed to wind 
stress relationship (discussed in Section 2), and these parameterization issues limit 
the accuracy of the stress fields. However, a greater limitation is posed by the long 
averaging periods, which result in a considerable loss of information on the temporal 
variability of the wind field. In addition to the wind stress fields derived from long- 
term ship reports, daily maps of wind speed derived from surface pressure analyses, 
supplemented by ship and buoy reports, are available. While these are available much 
more frequently than the averaged stress fields, they are not direct measurements of 
the surface wind field, but inferred from the much smoother quantity of pressure. 
Spatial resolution is considerably reduced due to these inherently large pressure field 
scales, and small, intense circulations may be lost in the broader scale unless they are 
observed independently. 

The insufficient amount of wind measurements, sparse global coverage, and re- 
stricted scales of conventional wind information available limit the studies that can 
be done. In particular, conventional sources generally provide winds only at larger 


scales (> 1000 km), geostrophically from surface pressure fields, and at very small 





scales (< 10 km), from ship and buoy reports. There is virtually no resolution of 
the intermediate spatial scales, which contain significant energy and are important to 
atmospheric forcing of the ocean (Freilich and Chelton, 1986). The conventional data 
also lack sufficient temporal resolution for modeling the temporal variability of the 
upper ocean, especially over larger areas. The absence of information at shorter time 
and space scales forces reliance on long time-averaged data for modeling purposes. 
These averages are of some use in studying the mean circulation, but lose much of the 
variability needed to drive the more complex models. Some studies of the temporal 
and spatial variability in realistic wind fields, as well as of the impact of that variation 
on ocean motions, have been conducted, though the studies are necessarily limited 
by the coarse resolution of climatological wind fields (for example: Willebrand, 1978; 
Willebrand et al., 1980; Muller and Frankignoul, 1981). The forcing effect of wind 
stress on the ocean surface is a major determinant for ocean waves and upper ocean 
currents. Knowledge of these motions is crucial to many fields, including basic oceano- 
graphic research, climate studies, and practical ocean operations. The interrelation 
of these motions and the wind forcing can be more fully evaluated only if an accu- 
rate, global wind field is available for input into the numerical models. Conventional 
sources of wind field information, including many different types of observations and 
analyses, do not currently provide sufficient information for this purpose, and a re- 
peated, global measurement of wind stress (or velocity) over the ocean would therefore 
be a unique contribution (O’Brien, 1982). 

The microwave scatterometer, yielding wind velocity as inferred from surface 
roughness measured by reflected radiation, was developed to provide this missing in- 
formation. Mounted on aircraft or, more effectively, on satellites, this instrument has 
the potential to provide a more accurate, variable, global wind field at higher temporal 
and spatial resolutions than currently available. The SEASAT Satellite Scatterometer 
(SASS), launched in 1978, was the first operational scatterometer system to be flown. 
The three months of data returned from that mission is a unique set of information 
that has yet to be repeated, since no scatterometer has flown since. Details of basic 


scatterometer theory, the SEASAT mission, and SASS itself are given in Section 2. 





The wind forcing information derived from scatterometer measurements provides 
both new research opportunities with the expanded information it makes available, as 
well as very useful information for immediate use in practical operations. Improved 
specification of the wind forcing would allow ocean modelers to better estimate ocean 
behavior, and therefore reach a more thorough understanding of upper ocean physics. 
The scatterometer can resolve winds on the intermediate spatial scales and the global, 
shorter time scales that are missing in conventional data. This information would con- 
tribute greatly to modeling of the temporal and spatial variability of the upper ocean 
(Freilich and Chelton, 1986). The ocean is forced by the wind stress, both directly, and 
indirectly through mass field adjustments, and the stress field must be incorporated 
into any realistic ocean model. Current parameterizations of thermodynamic forcing, 
including the surface fluxes of sensible and latent heat, use aerodynamic bulk for- 
mula models that incorporate wind velocity. So, better wind field information would 
improve our understanding and modeling capability of both wind stress and ther- 
modynamic forcing. In the absence of good forcing fields from conventional sources, 
oceanographers are restricted to artificial forcing fields, at best adjusting the sparse 
information currently available to obtain reasonable results from the model. Some 
of the new ocean modeling opportunities opened up by the more accurate, extensive 
wind fields provided by the scatterometer include: 1) improved ocean basin models, 
2) improved wind forcing related process studies, 3) improved prediction and moni- 
toring of ocean climate phenomena such as El Nino, and 4) better input for practical, 
daily ocean forecasting models. 

Scatterometer derived wind fields will also have valuable applications to meteo- 
rological research and weather prediction, and to a better understanding of processes 
at the air-sea interface. Assimilation of scatterometer data may have a significant 
impact on numerical weather prediction, through both model initialization and sub- 
sequent updates (Cane et al, 1981; Yu and McPherson, 1979; 1981; Duffy and At- 
las, 1986; Harlan and O’Brien, 1986). The additional information provided by the 
scatterometer is particularly important in data sparse areas, such as the Southern 


Ocean. Scatterometer derived wind fields also have important applications to naval 





and commercial ocean operations. These include tactical weather forecasting, acous- 
tic conditions forecasting, surf forecasting for amphibious operations, ship routing, 
and warnings of potentially catastophic situations. More commercial uses include oil 
platform design, drilling schedules, and commercial fishing schedules. The potential 
cost savings associated with these scatterometer wind applications are summarized in 
O’Brien (1982). 

Better vector wind stress fields will be needed to drive the increasingly sophis- 
ticated ocean models, designed to resolve temporal and spatial variability on finer 
scales. The choice of a particular wind field for this forcing is crucial when we want 
to more accurately establish the influence of the “real” surface winds on ocean cir- 
culation. At that point it becomes crucial to pick the “best” wind field available, 
in terms of accurately reproducing the actual winds with sufficient resolution at the 
desired scales. Even if it is not possible to determine which is best in an absolute 
sense, it is important to at least reach a good understanding of the differences be- 
tween the various wind fields available, especially in their ultimate impact on model 
output. In particular, when conventional and scatterometer derived fields are consid- 
ered as alternative forcing functions, the presence or absence of significant differences 
between the two can influence the decision to deploy a scatterometer in the first 
place. If no significant differences are found between the two fields, over the time 
and space scales of interest, then the effort and expense necessary to develop and 
maintain an operational scatterometer is not justified, since it would not significantly 
augment information already available. This study does not intend to address the 
more comprehensive question of the need for a scatterometer, assuming, as discussed 
briefly above, that this requirement is already well-established (Brown, 1983;1986; 
O’Brien, 1982). Rather, assuming that the scatterometer wind fields will be avail- 
able and will therefore compete with conventional winds as model input, this work 
seeks to examine the differences between these fields in a more detailed fashion. This 
is accomplished by qualitatively and quantitatively comparing two such fields, the 
conventional surface wind analyses from the Fleet Numerical Oceanographic Center 


(FNOC) in Monterey CA, and SEASAT Scatterometer winds from R. Atlas at the 





Goddard Space Flight Center (GSFC) (Atlas et al., 1987), obtained through the Jet 
Propulsion Laboratory (JPL), Pasadena, CA. An overview of satellite scatterometry 
is provided in Section 2 to introduce some basic atmospheric boundary layer and mi- 
crowave backscatter concepts. Section 3 describes the SEASAT mission and SASS, 
and Section 4 provides a general! description of the various data used and outlines the 
objective mapping procedure used to regrid the scatterometer winds. Section 5 pro- 
vides more detailed wind field statistics and various quantitative comparisons between 
the scatterometer and conventional winds, including boxed mean differences, scatter 
plots of various quantities, difference field statistics, and wind stress cur! calculations. 


Section 6 contains the summary and conclusions. 





2. Background 


The development of microwave scatterometer theory is based on studies of the 
atmospheric boundary layer, radar scatter, and ocean surface motions. Some of the 
basic concepts from these areas that apply directly to scatterometer work are reviewed 
here. The complex nature of the air-sea boundary makes quantifying and measuring 
physical processes there a difficult problem. Wind stress, 7, represents the transfer 
of momentum from the atmosphere to the ocean. This transfer of momentum slows 
the air near the ocean surface, creating an atmospheric boundary layer. The vertical 
profile of wind velocity in the atmosphere is dominated by molecular processes in 
the few millimeters closest to the sea surface, with turbulent motions, governed by 
friction and buoyancy, dominating the rest of the surface layer (Stewart, 1985). At 
higher levels the surface effect diminishes and the geostrophic general circulation of 
the atmosphere governs the flow. The region most critical to scatterometer work is the 
20 meters closest to the surface, and concentration on this region affects assumptions 
made about the relative dominance of various processes. In general, the wind velocity 
profiles are derived analytically, combining a logarithmic layer with an Ekman spiral 
modified by secondary flow for winds greater than 5 m/s (Brown, 1986). The various 
parameters and details behind this theory are discussed in more detail below and in 
Appendix 1. 

The study of atmospheric and surface boundary layers is based on formulas and 
empirical correlations that were originally derived from empirical studies. Many field 
experiments have been conducted to investigate the relevant processes, and then to 
adjust the formulas and refine the empirical constants (see Stewart, 1985 for partial 
review). Three quantities, the fluxes across the sea surface of horizontal momentum, 
sensible heat, and latent heat, are of particular interest in the 1 cm to 10 m surface 
layer where the surface effect dominates. Within this turbulent atmospheric surface 
layer, the fluxes are approximately constant with height and the wind and temperature 
profiles are logarithmic - with adjustments for stability (Halbertstam, 1980). Much 


study has been devoted to deriving the fluxes from wind speed and temperature 
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profiles. The similarity concepts of Obukhov (1946) and initial experiments of Monin 
and Obukhov (1954) are incorporated in most of the profile formulas. The vertical 
profile of wind velocity, and its relation to surface drag, or stress, strongly influences 
the transfer of energy from atmosphere to ocean. Wind stress, 7, is accepted as the 
major driving force for large-scale ocean circulation. However, practical considerations 
have dictated that, rather than being measured directly, the wind stress is normally 
derived from wind velocity, U, usually measured at some distance above the ocean 
surface. This is a highly complicated relationship, dependent on the complex dynamics 
of the boundary layer at the air-ocean interface. 

The scatterometer uses a measure of the surface roughness, determined from the 
reflection of microwave radiation, as an inferred measure of the air-sea energy transfer. 
The surface roughness can be connected to a particular geophysical parameter, such 
as wind velocity or wind stress, as well as to the backscattered power that forms the 
scatterometer signal. In practice, a direct relation is sought between the backscat- 
ter and a geophysical parameter, commonly chosen to be the wind velocity at some 
height above the surface. However, this relation has no direct physical basis, and is 
determined empirically. The theoretical justification for the empirical relation comes 
rather from our understanding of various parts of the process, and the interconnec- 
tion of quantities such as the wind velocity, wind stress, ocean wave spectrum, and 
backscattered radiation. The review of these basic concepts below is divided into 
three broad areas, the U-to-7 relation, the impact of wind stress on the ocean wave 
spectrum, and the relation between the ocean wave spectrum and reflected radiation. 
These are examined in turn, followed by a discussion of their practical, but indirect 
application in the direct relation between wind velocity and reflected radiation. In a 
sense this discussion can be seen as describing two paths to the same goal, the first 
a multi-step, theoretical process and the second a single, empirical step, where the 


latter path is dependent on the principles established by the former. 
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Wind velocity and wind stress 


The U-to-7 relation is determined by converting the wind at some reference 
height, U,es, through the boundary layer, to a surface stress, 7. This connection 
between the wind velocity profile and wind stress are of primary importance to scat- 
terometry. The interfacial momentum flux, or wind stress, can be defined as the 
average turbulent transfer of horizontal momentum from atmosphere to ocean by 
vertical air movements (O’Brien, 1982). For a constant wind, the wind stress is cal- 
culated from correlations between perturbations of the vertical and horizontal wind 


components, expressed as : 


where: 
7 = wind stress (N/m?) 
p = air density (kg/m°*) 
u,w = horizontal and vertical velocity (m/s) 


(Stewart, 1985) 


The averaging period must be long enough to form stable averages of the turbulent 
properties, yet short enough for mean conditions to be steady. Times of 40 minutes 
to 1 hour are usually chosen. Compounding the complex theoretical issues in this 
relation are the practical difficulties inherent in measuring the relevant quantities. 
The ocean surface is a difficult environment to work in, and the critical turbulent 
time and space scales are relatively short. To circumvent the constraints of direct 
measurement, a series of empirical relationships (bulk or aerodynamic formulae) have 
been developed. These relate the wind stress to the wind velocity at some height, 
with the usual form: 


T= pCyp laee 


where: 
Cp = drag coefficient 


U,e¢ = wind velocity at some reference height (m/s) 
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The exact formulation of the drag coefficient is still subject to considerable debate, 
and varies with wind speed and atmospheric stability, as discussed in Appendix 1. A 
commonly accepted formulation is that of Large and Pond (1981), which combines 
a constant value for Cp at lower wind speeds (< 11 m/s) with an expression that 
is linearly dependent on wind speed at higher speeds (see Appendix 1 for the exact 
form). 

The wind velocity profile near the surface varies sharply with height and depends 
on both 7 and atmospheric stability. The transfer of heat and momentum depend on 
turbulence generated by both wind shear and buoyancy (affected by temperature and 
humidity stratification). Shear dominates near the surface and buoyancy becomes 
relatively more important at some height above the surface. The influence of stability 
is illustrated by considering the case of cold air blowing over warmer water. Air is 
heated from below and rises, causing instability in the air column and enhancing 
the turbulence due to just the mechanical shear of the wind stress. In the opposite 
case, when warm air blows over cold water, there is a suppression of turbulence 
as cooling from below stably stratifies the air column. The impact of wind stress 
and stability on motions in the turbulent layer is described quantitatively by Monin- 
Obukhov similarity theory, discussed in more detail in Appendix 1. A primary result 
from these theories is that, assuming that the influence of atmospheric stability is 
weak, the wind velocity profile is logarithmic with height for neutral stability, and 
simply shifted by stability corrections. Under unstable conditions, or for a strongly 
stratified atmosphere, these formulations are no longer valid, but the condition of 
weak stratification is generally satisfied over widespread areas, except in the tropics 
where strong convection cells are common. 

Though considerable debate still exists on some details of the theory and on the 
exact formulation of the relation, there is a generally accepted understanding of the 
wind stress-to-velocity connection. The choice of wind velocity in particular, and if 
so, at what height, as the relevant parameter to relate to microwave backscatter is 
a still greatly debated issue. The neutral stability wind at 19.5 m was chosen for 


the original processing of SEASAT scatterometer data, based on anemometer heights 


13 





aboard some weather ships. From boundary layer considerations, a better choice 
would be the wind velocity closer to the surface (such as at 10 m), or the friction 
velocity. Another alternative parameter, using the wind just above the surface, is 
discussed in a later section. Though the “best” parameter has yet to be decisively 
chosen, for the initial evaluations made for SASS development, the most practical 
choice was wind velocity. The boundary layer wind velocity theories described here 
provide justification for that particular choice, as they connect the wind velocity 
to wind stress. A factor requiring further evaluation is the effect of stratification, 
especially since it can have a significant influence at all scales. The parameterizations 
discussed here hold for relatively weak stratification effects, but are questionable under 
conditions where atmospheric stratification plays an important role, such as in the 
tropics. To isolate these effects for future uses, widespread measurements are needed 
of both the currently available sea surface temperature measurements, and the not yet 
generally available air temperature measurements over the ocean (Atlas et al., 1986). 
A more detailed parameterization of the effects of stability, along with better sea 
surface and air temperature observations would significantly improve and extend the 


range of scatterometer wind algorithms. 
Wind stress and the ocean wave spectrum 


Given some stress at the ocean surface, the next step is to determine how this 
transfer of momentum from the air is reflected in the ocean wave spectrum. Un- 
derstanding the processes behind this energy transfer is a fundamental problem of 
remote sensing. Important factors in the wind-wave relation include wave generation, 
Wave-wave and wave-current interactions, wind strength, atmospheric stratification, 
sea surface temperature, and sea state. The ultimate goals are to determine how 
much energy is transferred from air to ocean, on what scales the transfer takes place, 
and how the energy is partitioned to various ocean features, in particular, waves and 
currents. Wind generated waves range from millimeters to hundreds of meters and, 


while there is considerable theoretical understanding of the longer gravity waves, the 
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high-frequency end of the spectrum is much less well understood. Current theory 
holds that the short wavelengths reflect the surface response to the local wind field. 
As discussed in detail in the following section, scatterometers are designed to look 
at these wavelengths, in particular at the centimeter-length capillary waves. Some 
physical factors that influence the surface ripple field are: 1) the tangential and nor- 
mal stresses at the surface, 2) the fraction of momentum transfer going into capillary 
waves, 3) the magnitude and direction of the swell, 4) ocean currents and mixing, 
5) wave history, 6) hydrostatic stability (Davidson et al., 1981). The lack of definitive 
descriptions of open ocean capillary waves and their relation to the surface fluxes 
has hindered theoretical understanding in this area. Evidence also exists that only a 
portion of the air-sea momentum transfer goes into these waves, and the assumption 
that the short waves reflect the entire momentum transfer, and therefore the wind 
stress, is not supported by any complete quantitative theory (Brown, 1986). Partial 
explanations are available for various processes, especially wave generation, growth, 
and interaction, but many critical effects have yet to be modeled, and a single, unified 
theory does not exist. Short waves are associated with skin friction-induced small- 
scale stress, and therefore do not account for the direct momentum transfer through 
form drag on long waves from pressure fluctuations (Brown, 1986). The partitioning 
of momentum transfer between shear stress and form drag is discussed in more de- 
tail in Appendix 2. Other unaccounted for factors that may distort the local wave 
field to local wind relationship include swell, usually generated at more distant loca- 
tions and unrelated to local winds, and momentum transfer to ocean currents, thus 
unavailable for wave generation. At present, the influence of uncorrelated swell and 
the momentum loss to currents are assumed to be small and are neglected. A short 
review of some theories and experiments on the basic physics behind the transfer of 
energy from air to water is given in Appendix 2. The concepts and studies discussed 
here provide the still rough theoretical understanding of the wind-wave relation used 
to justify scatterometer development. These present interpretations are supported by 
the generally good results gained from scatterometer measurements made thus far. 


Further refinements of the assumed stress-to-wave spectrum relation are still needed 
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to resolve effects such as the time response of wave formation, interference from other 


waves, and the effect of sea slicks. 
Ocean wave spectrum and radar backscatter 


The final step is to determine the ocean wave spectrum from information on 
surface roughness gained from reflected radiation, or backscatter. The interaction of 
electromagnetic radiation with the surface provides information on surface features, 
depending on the wavelengths of the radiation and the features. Extracting surface 
information from the power of the backscatter signal therefore requires knowledge of 
both scattering theory and surface characteristics. Details of relevant surface scat- 
tering principles are provided in Appendix 3. The primary radiation measure is the 
normalized radar cross section (NRCS), o° , a ratio of the reflected to incident energy 
across a unit surface area. The crucial problem is to characterize the shape of o° 
for a given sea surface parameter, generally chosen to be the ocean wave spectrum. 
Scatter from the wavy ocean surface is described by two physical mechanisms, spec- 
ular and Bragg (or resonant) scatter, which depend on the incidence angle (angle 
between the radar beam and local surface vertical) and wavelength of the incident 
radiation, and on the wavelengths of the surface waves. Appendix 3 contains further 
descriptions of these scatter mechanisms, including the incidence angle dependence 
and ranges, and the adjustments needed to accomodate the multiple scales of oceanic 
motion. The most useful return for scatterometry is the resonant scatter from capil- 
lary waves, yielding both speed and direction information on the wave spectrum, and 
by association, on the surface stress or wind velocity. 

These basic scattering concepts are complicated by the complex conditions actu- 
ally present on the sea surface. Wave tank and field experiments have shown spikes 
in backscatter returns due to sharp-crested waves that are close to breaking (Kwoh 
and Lake, 1984). The effect on o° of rain striking the surface is significant, but 
not well understood. The capillary waves tend to be damped by rainfall, thus de- 


creasing backscatter return, but rain-generated ripples also appear, increasing the 
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radar return. Some experimental evidence has shown go” to increase with rainfall 
(Moore et al., 1979), but more study is needed to determine the magnitude of this 
effect. Sea slicks can significantly affect backscatter as well (Huhnerfuss et al., 1983), 
and recent work further demonstrates the impact of wave slope and atmospheric strat- 
ification on the radar return (Keller et al., 1985). At higher wind speeds, when the 
surface is confused, with significant foam and spray coverage, the application of Bragg 
scatter theory is questionable and new models may be necessary {Atlas et al., 1986). 

While microwave backscatter theory is reasonably well understood in terms of 
electromagnetic radiation principles, the precise connection to the short wave spectra 
is not well specified. Rather, as discussed in detail in Section 3, the assumed relation 
between the two is used to justify empirical models that relate backscatter to wind 
velocity and various parameters of the radar signal. Though these models produce 
realistic results and have been used with much success, they shed little light on the 
physical processes that link microwave backscatter and ocean wave spectra. Also, as 
discussed by Plant (1986), the current formulations of short wave spectra derived from 
physical principles (eg: Phillips, 1977; Kitaigorodskii, 1983; and Phillips, 1985) are 
not directly applicable to scatterometry, since they cannot account for the dramatic 
increase in the backscatter signal with the increase in wind speed (and therefore sur- 
face roughness). Addressing this question, Plant (1986) derives a form for the short 
wave spectrum that incorporates long wave-short wave interactions as well as wind 
input and dissipation. The expression for the short wave spectra is valid to second 
order in long wave slope, and is then used, for the case of locally generated seas, to 
develop an expression for 0? , also valid to second order in long wave slope. This 
formulation allows a more detailed investigation of the effects of long wave slope and 
propagation on backscatter, through the tilting and hydrodynamic modulation of the 
short waves, but it is still a considerable simplification of actual ocean conditions. 
Assuming that the relevant parameter is the backscatter from capillary waves, in the 
real ocean these ripples (and therefore the associated o° signal) are modulated by 


ocean currents,internal waves, and sea state, and subject to tilt effects from longer 
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waves. Surface currents cause variations in o° through their effects on waves, includ- 
ing changing the direction of propagation and amplifying wave amplitudes by local 
velocity convergences. The effect of internal waves is not that pronounced at this 
scale, but sea state can have a significant impact (Atlas et al., 1986). The influence of 
sea state is especially important in rougher conditions, when foam, spray, and break- 
ing waves disrupt the surface and contaminate the backscatter return, often as signal 
spikes. 

In the absence of detailed knowledge of the actual long wave field, the optimum 
conditions for scatterometry occur when equilibrium is reached by those waves approx- 
imately one order of magnitude longer than the capillary waves (Atlas et al., 1986). 
These conditions are assumed in scatterometer algorithm development, and though 
this steady state condition is typical for the open ocean, it will not hold near signifi- 
cant circulation features, such as fronts. In order to completely specify the modulation 
and tilt effects for future algorithms, more extensive and more highly resolved itn situ 
ocean wave spectra information is needed. In addition to refining the empirical scat- 
terometer relations, this information could also provide a means of developing and 
validating theories on the physical processes discussed here. These scattering con- 
cepts, though simplifications of realistic conditions, are supported by the few exper- 
imental studies available (see Stewart, 1985 for summary). The difficulty of making 
simultaneous measurements of ocean wave spectra and backscatter has necessarily 
limited the amount of experimental evidence. While many questions still remain to 
be answered, there is increasing physical understanding of microwave backscatter, of 
the short ocean wave spectrum, and of the connection between them. Future study 
will be directed to quantifying the effects on microwave backscatter of rain, wave 


slope, sea surface temperature, currents, internal waves, and breaking waves. 
Empirical U-to-o° relation 


The theories described above are used to establish a physical foundation, through 


a convoluted, multiple-relation process, on which to connect microwave backscatter 
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surface wind velocity. ‘he response of short capillary waves to environmental changes 
at the air-sea interface, as measured by radar backscatter, is used to determine wind 
speed and direction. Given the theoretical support provided by these concepts, the 
practical application is to use that basis to justify seeking a direct relation between 
wind and the measured radar cross-section. The direct relation is established in three 
steps: 1) isolate the critical parameters, 2) establish the basic form of the equation 
(ie: power law, drag law, etc.), and 3) describe the detailed aspects of the model 
through empirical coefficients. Empirical studies have shown the dependence of o° 
on wind speed and wind direction relative to the horizontal antenna pointing angle 
as well as on the incidence angle, polarization, and frequency of the radar signal 
(Plant, 1986). Wind speed dependence is described by an adjusted power law while 
directional dependence is characterized by an approximately cos(2x) dependence on 
the relative azimuth angle, x. The empirical coefficients are derived by tuning the 
basic equation to a set of data, for which both the known and unknown variables have 
been found. An obvious limitation to this process is that tuning to a particular data set 
constrains the backscatter-to-wind relation to the range of conditions present in that 
data. The specific details of the 0° -to-wind relation used for SEASAT scatterometer 
data are described in Section 3. It should be mentioned that the radar cross-section 
has been shown to correlate well with both wind velocity and friction velocity (Liu 
and Large, 1981). However, as mentioned above, due to the need for conventional 
measurements to validate the empirical relation, and to the relative unavailability of 
wind stress measurements, a wind velocity dependence was chosen for SEASAT. One 
problem with current formulations is an inconsistency between the parameterizations 
for neutral wind and friction velocity. The current choices of power law relations to 
backscatter for both are not compatible with the expected connection between the 
two for a neutrally stratified atmosphere (Pierson et al., 1986). 

Recent theoretical and experimental evidence indicates that the capillary wave 


spectrum, and by assumption the microwave backscatter, correlates better with a 
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near-surface wind field characteristic, R, than with either the neutral wind or friction 


velocity (Pierson et al., 1986). This new parameter is expressed as: 


_ Be 





where: 
U(4) = average wind at half a Bragg wavelength above the 
surface 
C(A) = phase speed of Bragg waves 
A = Bragg wavelength 
(Pierson et al.,1986) 


This term is related to the normal stresses, and correlation with short wave spectra 
depends on the assumption that normal stresses dominate tangential stresses in wave 
generation (Atlas et al., 1986). Unlike a direct correlation of backscatter to friction 
velocity, use of R requires definition of the wind profile and specification of a drag 
coefficient, or friction velocity and roughness length (Pierson, et al., 1986). So the ap- 
plication of R with scatterometer data is similar to using the neutral stability winds, 
but with a characteristic that is more physically related to the surface roughness. 
These concepts may improve future scatterometer missions, but SEASAT data was 
processed with a o®% -to-U relation and that will be emphasized here. Despite many 
unresolved issues, the available SASS results have demonstrated that useful correla- 
tions between microwave backscatter and a wind or surface roughness parameter can 


be found (Brown, 1986). 
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3. SEASAT Scatterometer (SASS) 


A considerable microwave backscatter history culminated in the development 
of the SEASAT scatterometer, commencing with experiments on radar sea clutter 
during and after World War II. The initial proposal to use microwave backscatter 
to determine winds at the ocean surface came from Moore and Pierson (1967). A 
history of the various initial experiments can be found in Moore and Fung (1979). 
Following the general procedures discussed above, various empirical forms were devel- 
oped for the backscatter-to-wind vector relation (Moore and Fung, 1979; Boggs, 1981; 
Jones et al., 1982; Schroeder et al., 1982a; Woiceshyn et al., 1986). The final choice 
for SEASAT, discussed in detail below, was a single power law relation with a least 
squares inversion (Moore and Fung, 1979; Jones et al., 1977). Following the SEASAT 
mission, the scatterometer model function was tuned through comparison with high 
quality surface wind data, primarily that collected in the Gulf of Alaska (GOASEX) 
and Joint Air-Sea Interaction (JASIN) experiments (Jones et al., 1982; Schroeder 
et al., 1982a, 1982b; Wurtele et al., 1982; Brown, 1983). Some results from these 
experiments are summarized following a more detailed description of the SEASAT 
scatterometer. 

The SEASAT oceanographic satellite, launched in June 1978, carried several 
microwave remote sensing instruments, including SASS (Figure 1). The mission of 
SEASAT was to conduct proof-of-concept experiments for the detection of surface 
features by remote sensing. These features included surface wind, sea surface temper- 
ature, wave height, ocean topography, and sea ice, as well as information on internal 
waves, atmospheric water vapor, and the marine geoid. The SEASAT mission ended 
in October 1978 due to a power failure, but the three months of returned data provided 
sufficient information to meet the original objectives. The satellite flew in an approx- 
imately circular orbit with an inclination angle of 108°, at an altitude of 800 km, and 
with an orbital period of about 101 minutes, thus circling the Earth 14 times a day 
(Boggs, 1982). The wide swath scatterometer sensor covered 95% of the ocean every 


36 hours (Boggs, 1982). During the first part of its mission SEASAT maintained 
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a 17-day repeat cycle (the period before repeating the same ground track), with a 
minimum spacing between equatorial crossings of 88 km during that interval. From 
August 27, 1978 to the end of the mission, SEASAT switched to a three day repeat 
cycle, with a minimum equatorial crossing separation of 470 km (Boggs, 1982). The 
scatterometer operated nearly continously during the mission, returning over 90 days 
of good information. SASS used four, dual-polarized fan beam antennas, oriented 
at 45° and 135° relative to the satellite track, to get measurements at two orthog- 
onal azimuth angles. The X-shaped illumination pattern of these beams formed a 
double-sided swath of potential wind-vector data about the satellite subtrack, and a 
narrow strip of only wind speed information near nadir (Figure 2). The transmitted 
frequency was 14.6 GHz, with a slight shift in the return signal from doppler effects. 
The antenna beams were electronically subdivided into 15 resolution cells by doppler 
filtering, using the intersection of the antenna pattern and the lines of constant doppler 
shift. The final backscatter resolution cell was generated by integrating the received 
power over a 1.89 second period. The signal processing provided a spatial resolution 
of approximately 50 km within the double-sided swaths. 

Several sources of instrument error affect the backscatter measurement, 0° , in- 
cluding random errors due to communication noise, uncertainty in spacecraft attitude, 
and instrument processing, as well as several bias errors (Boggs, 1982). Communica- 
tion noise is the primary random error and depends on the doppler bandwidth, the 
signal integration period, and the signal-to-noise ratio (Fischer, 1972). The next most 
significant source of measurement error comes from uncertainty in spacecraft orienta- 
tion, and therefore antenna pointing direction. The attitude error encompasses errors 
in roll, pitch, and yaw, as well as instrument alignment. Boggs (1982) lists values for 
these various random error sources. Bias errors between the four antennas were esti- 
mated from measurements of backscatter over the Amazon rain forest and relatively 
small corrections were made during sensor data processing. Corrections to o° for at- 
tenuation through the atmosphere were made for backscatter cells that coincided with 
measurements made by the Scanning Multichannel Microwave Radiometer (SMMR). 


This was possible only on the right side of the satellite track, since the SMMR swath 
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was one-sided. The lack of attenuation correction is significant primarily for high 
rain rates, and the anomalous wind values produced under those conditions are easily 
flagged in most cases. Finally, the 0° measurements from the integrated “footprints” 
were screened for anomalous values and for the presence of ice or land, and then input 


into a geophysical algorithm to obtain the wind vectors. 
SASS geophysical algorithm 


The particular geophysical algorithm used for SASS is described completely in 
Boggs (1982). It computes the 19.5 m neutral stability wind speed and direction 
at a particular time and location. Basically, it consists of three components: 1) 
a cell-pairing process to match forward and aft cell values, 2) a wind-to-o° model 
function, and 3) a least-squares estimator to invert the model function. Prior to the 
backscatter-to-wind inversion, a complex temporal and spatial reorganization of the 
backscatter data is needed. The o° cell measurements are spread out across the swath 
for all particular times and antenna patterns. The wind determination algorithm 
requires at least two, roughly colocated (in both time and space) o° measurements, 
approximately orthogonal in azimuth. The orthogonality condition requires at least 
one o° measurement from both the forward and aft antenna beams. The data can be 
grouped either by cell-pairing, in which individual forward and aft cell measurements 
are matched, or binning, in which all o° measurements within a given area are used 
to find a single wind solution for that area. Cell pairing was used in the production 
of the original wind vector data produced from SASS and is described more fully 
in the next section. Forward and aft beam measurements were paired if they fell 
within given time and distance separations, with some redundancy of data as some 
individual cells were paired more than once. The cell-pairing mode was chosen for the 
initial data production since it yielded the highest resolution wind vector solutions 
(Boggs, 1982). Subsequent processing (discussed at the end of this section) used the 
binning mode, with a 50 km resolution. In both cases, nadir wind solutions (in the 


narrow strip directly below the satellite) were formed from binned data, without the 
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need for orthogonal measurements, due to the lack of directional information at low 
incidence angles. In each case, the time and location assigned to the wind solution 
are the centroids of those of the grouped 0° measurements. Once matched, the paired 
or binned backscatter measurements are input into the model function, which is the 
empirical relationship that describes the dependence of o° on the 19.5 m neutral 
stability wind. The relation chosen for SASS data can be expressed in logarithmic 


form as: 


Cap = G(4, x,€) ae H(8, x, €)logioU19.5 


where: 
oS, = backscattered power in decibels 
@ = incidence angle 
x = wind direction relative to radar azimuth angle 
€ = polarization of incident radiation 
Ui9.5 = neutral stability wind at 19.5m (m/s) 
(Boggs, 1982) 


The o° values have a range of more than 50 db for typical wind speed ranges, 
so there is a clear advantage to using a logarithmic form (Woiceshyn et al., 1986). 
It should also be repeated that this particular empirical relationship is an arbitrary 
choice and does not have a direct physical basis. The incidence angle, 6, is defined 
at the surface as the acute angle between the antenna look direction and the local 
surface vertical at a doppler cell center (Figure 3). The relative wind azimuth, y, is the 
angle between the local horizontal wind direction, -y, and the radar azimuth angle, 
@, measured from north to the projected antenna look direction, at the cell center 
(Figure 4). The radar azimuth angle can be thought of as the instantaneous azimuthal 
direction of the antenna beam pattern on the surface (Woiceshyn et al., 1986). So, the 
model function depends implicitly on the actual wind direction through an explicit 
relation to the relative wind azimuth angle, y. y is defined so that the upwind 
condition, with the antenna pointing directly into the wind, corresponds to xy = 


O° and the downwind condition to x = 180°. The original sensor data processing 
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actually used the value of the radar azimuth angle calculated at the sub-satellite 
point, instead of at the cell center. This introduces a relatively small error that 
can be neglected in most cases (Boggs, 1982). The wind direction used for SASS 
processing is defined in the meteorological sense as coming from a given direction, 
measured clockwise from the north. The empirically derived G and H coefficients are 
given as coefficient tables for vertical and horizontal polarizations and incremental 
values of 6 and y. They are tabulated for 0° < @ < 70° in 2° steps and 0° < x < 180° 
in 10° steps for both polarizations. A first-order interpolation was used to extract 
values from the 2° @ table, while a second-order interpolation was used with the 10° x 
table. The dependence of G and H on azimuth is approximately cos(2x*), where x* 
is the difference between wind direction and radar azimuth, with maxima at y* = 0° 
and 180°. The harmonic dependence means that the relative wind direction, y, used 
in the model function, has a functional range of 180 degrees (Boggs, 1982). The G-H 
tables were originally determined from aircraft measurements and then modified after 
several SEASAT validation workshops. The final table was formed by combining two 
models over different incidence angle ranges (Boggs, 1981; Schroeder et al.,1982a). 
The selection of parameters and the choice of a relatively simple, logarithmic form 
with a power law dependence are still subject to considerable debate (Atlas et al., 1986; 
Woiceshyn et al., 1986), but generally good results and confirmation of basic concepts 
were obtained from this initial formulation. The model function is a unique relation 
for the dependence of o? on the 19.5 m neutral stability wind, but the wind retrieval 
process requires using this relation in the opposite direction. This inversion requires 
at least two colocated, roughly orthogonally viewed o° measurements and is not a 
unique specification. A weighted least-squares estimation, also referred to as a sum- 
of-squares approach, was used to accomplish the actual inversion (Boggs, 1982; Jones 
et al.,1982). 

One limitation of this o° -to-wind algorithm is that for a given o° grouping 
there may be up to four possible wind vector solutions, all with approximately the 
same wind speed, but with widely varying directions. Figure 5 illustrates this for 


the ideal case of two orthogonal, noise-free, colocated 0° measurements of the same 
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wind. The curves are possible solutions in wind speed (U)-azimuth angle (y) space 
for particular values of backscatter (o° ) and incidence angle (@;). Wind solutions are 
found at the curve intersections and this case shows the large variation of directions 
and small range of wind speeds for the different possible solutions. As shown, a 90° 
separation in azimuth yields the most distinct intersections. Noisy measurements 
further complicate this process, since the intersections are then spread over a larger 
area. These direction aliases cannot be removed without further processing, usually 
requiring some independent knowledge of the wind field. Future scatterometers will 
help resolve this ambiguity problem, by adding an additional antenna on each side, 


reducing the maximum number of solutions to only two, separated by nearly 180°. 
SASS performance evaluations 


Evaluations of the final SASS wind retrieval algorithms initially indicated that it 
met the required specifications of +2 m/s in wind speed and +20° in wind direction 
on a synoptic scale. The final results from the SASS workshops were rms differences 
of 1.3 m/s for a range of 4 to 26 m/s, and 16° for a 0°-to-360° range (Lame and Born, 
1982; Jones et al., 1982). However, these rms values were calculated for a restricted 
data set that was limited to the particular wind speed range and geographical area 
of the conventional wind data used for tuning and validation. These statistics do not 
always accurately reflect performance for other subsets of the data, as discussed in 
detail by Woiceshyn et al. (1986). Sources of error for the SASS generated winds 
exist both in the geophysical interpretation of the radar signal and in factors that 
may contaminate the signal, such as rain, sea surface temperature, and radar polar- 
ization. Many of these questions were addressed in the scattering theory discussion. 
An additional point is that the inherent spatial averaging of each SASS footprint may 
smooth out variations of the small scale turbulent winds (up to 10 m/s) embedded 
in the footprint area (Brown, 1986). The spatial smoothing may be especially signif- 
icant for small, intense features, such as thunderstorms. Strong circulation features, 


such as fronts, cyclones, and hurricanes, also pose a problem, since they generally 
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violate the steady state conditions assumed in the o° -to-U algorithms. Ditferences 
in wind retrieval system behavior across a front are evident in the changed charac- 
ter of the scatterometer signal across the front (Brown, 1983). Adjustments to the 
wind retrieval system are needed to more fully account for these flow features, most 
likely as parameterizations that can be incorporated into the steady state algorithms 
(Brown, 1986). 

Various evaluations of the SASS data have identified problems with different 
aspects of the wind retrieval algorithms. Woiceshyn et al. (1986) use intercomparisons 
of horizontally and vertically polarized SASS data and comparisons between itn situ 
and SASS winds to demonstrate the inadequacy of the power law relation at all 
observed wind speeds and incidence angles. This stems from the tuning of the data 
to the limited range of conditions observed in the “surface truth” experiments and to 
the inability of a single power law to model the difference in surface roughness over 
two distinct wind regimes. Woiceshyn et al. (1986) discuss other problems, including 
the inaccuracy of the SASS winds at low wind speeds. The combination of a sum-of- 
Squares inversion technique and the signal from noise subtraction process may yield 
negative backscatter values for light winds (Pierson et al., 1986). Using a power 
law assumption with this inversion process means that the negative values must be 
discarded, causing a considerable loss of information at lower wind speeds. A better 
option may be to use a maximum likelihood estimation inversion that will make use 
of all the data (Pierson et al., 1986). Preliminary work on a new scatterometer wind 
extraction system that removes many of these biases in the SASS model is presented 
in Woiceshyn et al. (1984). A careful statistical analysis of the original SASS winds 
by Wentz et al. (1984, 1986) also revealed several systematic errors in the SASS data. 
They found that horizontally polarized backscatter signals yielded winds that were 
biased high relative to those from vertically polarized backscatter returns, a cross- 
swath wind gradient error was introduced from an incorrect 0° vs incidence angle 
relation, low signal-to-noise ratio o° values were discarded - so low winds were biased 
high and data gaps existed for areas of light winds, and the SASS winds were biased 


high by approximately 1 m/s relative to conventional winds. Some of these errors are 
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similar to those discussed by Woiceshyn et al. (1986), but Wentz (1986) retains the 
power law relation and corrects the model coefficients, rather than choosing another 
form for the model function. Finally, the question of directional ambiguity in the wind 
vector solutions has received considerable attention (Wurtele et al., 1982; Peteherych 
et al., 1984). Several alternatives to the SASS algorithm method have been developed 
(Hofman, 1982, 1984; Gohil and Pandey, 1985; Woiceshyn et al., 1986; Brown, 1986), 


and the design of future scatterometers will considerably reduce this alias problem. 
SASS data reprocessing 


The original SASS wind data, processed with the model function and cell-pairing 
method described above, is available on Geophysical Data Record (GDR) tapes lo- 
cated at the Jet Propulsion Laboratory (JPL). The alias problem and other errors 
discovered in the initial GDR processing, mentioned briefly above, motivated a re- 
processing of the entire sensor data set. Wentz et al. (1986) reprocessed the entire 
96-day data set, binning the backscatter values into 50 km cells oriented perpendicu- 
lar to the satellite subtrack. These combined 0° measurements were then inverted to 
obtain the wind vectors, on both 50 and 100 km grids. The winds were retrieved with 
an improved 0% -to-wind model, based on an assumed Rayleigh wind distribution, 
rather than tuning to tn sttu measurements. The model was designed to minimize the 
systematic errors dependent on polarization or incidence angle identified earlier, and 
used a single power law relationship, except for an adjustment for the higher nadir 
winds. Based on climatology, the winds were assumed to be Rayleigh distributed 
about a mean of 7.4 m/s (Wentz et al., 1984). This statistical approach has several 
advantages over a tuning approach, including the ability to use all the satellite mea- 
surements, rather than just those matched with a limited, high-quality conventional 
data set, and elimination of the problem of matching the different temporal and spa- 
tial scales of satellite and conventional data. The pre-averaging of the backscatter 
measurements across larger cells increases the signal-to-noise ratio (SNR) for the final 


o° value, which is especially useful at low wind speeds. A typical cell in the middle 
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of the swath contained four 0? measurements from each antenna, and SNR increased 
by a factor of two (Wentz et al., 1986). These reprocessed winds show more con- 
sistency and lower residual errors than the original SASS winds in intercomparisons 
over different polarizations and incidence angles. A comparison to buoy winds showed 
good general agreement in wind speed, with a 1.6 m/s rms difference and -0.1 m/s 
bias, supporting the statistical assumptions made in deriving the new model (Wentz 
et al., 1986). 

The directional ambiguity in the SASS winds was removed by several different 
methods, both subjective and objective, for various subsets of the original data. Pete- 
herych et al., (1984) developed a subjective dealiasing method and produced a global, 
15-day (6-20 September 1978) set of dealiased SASS winds. Their method is based 
on procedures described by Wurtele et al., (1982) and Baker et al., (1984), where 
meteorological analysis and pattern recognition are used to select an alias. Rather 
than using subjective meteorological analysis to dealias the data, alternative objec- 
tive techniques have been developed, including those of Hoffman (1982, 1984), Baker 
et al., (1984), and Yu and McPherson (1984). Atlas et al., (1987), used a modified 
version of the model of Baker et al., (1984) to objectively dealias the complete 96-day 
set of reprocessed, 100 km resolution, multiple-solution vectors. Atlas et al., (1987) 
employed the Goddard Laboratory for Atmospheres (GLA) analysis/forecast system 
in an iterative, three-pass procedure, combining conventional data and previously 
analyzed SASS winds as a first-guess field. These unambiguous SASS wind vectors 
are subsequently identified as Atlas winds in this paper, to distinguish them from 
the original GDR winds. They form the principal set of scatterometer winds to be 


evaluated and compared with conventional wind data in this study. 
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4. Data description and handling 


The initial data set obtained was a subset of the GDR data, supplied by NASA 
Ocean Data Systems (NODS) at JPL. Figure 6 shows a swath of these winds from 
SEASAT revolution 184. The alias problem discussed above limited its usefulness, 
and when the Atlas data set became available it was used, rather than the GDR 
data, for all subsequent evaluations. As mentioned in Section 3, the Atlas SASS 
winds are referred to as Atlas winds, in order to distinguish them from the original 
GDR data and from other sets of SASS wind data. In the preliminary analysis, we 
removed the directional ambiguity for selected GDR winds subjectively, by comparison 
of the plotted GDR wind vectors to synoptic weather and pressure maps. The SASS 
solution in closest agreement with the conventionally derived wind field was selected 
wherever possible, while those cases where an unambigous solution was not obvious 
were dropped from further consideration. Due to uncertainties imposed by this rough 
dealiasing scheme, the GDR winds were most useful in an initial screening of the Atlas 
winds, rather than as input for further analysis. After basic pattern agreement was 
established, these winds were most useful in ensuring the proper reading, transfer, and 
display of Atlas wind vector data. This study concentrates on a one-month subset 
of the full 96-day record, from 12 August to 10 September 1978, over a portion of 
the western North Atlantic, from 20-30°N and 40-80°W. The Atlas data is regularly 
spaced along the satellite subtrack, but it is irregularly spaced on a latitude-longitude 
grid, and the nature of satellite orbit ground coverage means that data samples in 
a given area are widely separated in both time and space. Typically, two passes 
cross the study area, separated by approximately 100 minutes and 2500 km, followed 
12 hours later by another two passes, offset in space from the first pair. 

Since the primary purpose was to compare scatterometer winds to a conventional 
wind field, a readily available, suitable conventional wind field was needed. The 6-hour 
synoptic surface wind analyses provided by the Fleet Numerical Oceanographic Center 
(FNOC) in Monterey, CA, through the National Center for Atmospheric Research 


(NCAR), were selected. These winds are generated by combining ship and buoy wind 
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reports with pressure-derived winds and are inherently smoothed to some extent on 
larger pressure field scales. A subset of the analyses was extracted, covering the same 
western North Atlantic area as the scatterometer data, interpolated to a 1° latitude by 
1° longitude grid. In addition to this wind field data, synoptic weather and pressure 
charts were also employed. These small-scale, large-area facsimile charts are quite 
crude, but provide a source of qualitative information for comparison to both the 
Atlas and FNOC fields. 

An important factor in any comparison between satellite and conventional mea- 
surements is the different temporal and spatial averaging inherent in each. The scat- 
terometer footprint is at the shorter end of the mesoscale range, a measurement av- 
eraged over roughly 50 km in space and nearly instantaneous in time, but sampled at 
widely separated intervals of both space and time. Ship and buoy wind measurements 
are point measurements in space, and may be averaged over an arbitrary time period, 
especially in the case of the buoy winds. Conventional, pressure-derived surface wind 
fields have larger, synoptic scales, comparable to those of the scatterometer swath 
(500 to 1000 km). These differing scales are important to a proper interpretation of 
the differences between the various data (Pierson, 1983; Stewart, 1985). 

So the available data are the GDR winds, the Atlas winds, the FNOC winds, 
and the synoptic weather and pressure maps. The primary comparison is between 
the Atlas and FNOC winds, both between vector fields and in a broader, statistical 
sense, with the weather maps used mostly for initial screening and evaluation. This 
screening included verifying the general flow patterns and the location of significant 
meteorological features. The Atlas-F NOC wind comparison was hindered immediately 
by the irregular spacing of the Atlas data (relative to a latitude-longitude grid). Boxed 
mean differences were calculated for some areas, by calculating component means for 
each field for a series of areas, then finding the differences between the means in these 
boxes. These have limited usefulness, especially since the resolution of the Atlas- 
FNOC difference field generated in this manner is very poor, so further processing of 


the scatterometer winds was necessary to regrid the field. To obtain regularly gridded 
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scatterometer data, objective mapping techniques were employed, concentrating on 


spatial mapping. 
Objective mapping procedure 


Objective mapping is basically an interpolation scheme that incorporates knowl- 
edge of the field statistics and data locations in construction of an informed estimate 
of the field at desired locations. A major advantage of the process is that it provides 
an estimate of the error field at the same time, yielding an objective measure of how 
well the mapping does at different points. The formulation used here follows that of 
Bretherton, Davis, and Fandry (1976), and most details of the actual mapping proce- 
dure are omitted, except for the form of the data covariance matrix, the general form 
and particular choice of which are described briefly below. Input of field information 
is made through choices of the form and scales of the data covariance matrix. This 


matrix, A,,, has the general form: 


where: 
Ars = covariance matrix of all pairs of observations, r and s 
F(x, — x,) = covariance as a function of spatial separation 


E = variance due to measurement noise 


The covariance function, F(x), where x is the spatial separation in x and y, is used to 
form the covariance matrix between of pairs of individual data measurements. After 
experimenting with various forms of this covariance function, a simple gaussian form 


was chosen: 
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where: 
F(x) = covariance function 
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Rs = signal scale (for both x and y in this case) 


ox = “noise” variance due to small scale signal 
Rn = “noise” scale 


Both of these matrices, A,;, and F(x), are dependent only on the position of the 
observations, so they are constant for different realizations of the data field, unless 
one expects the statistics of the data field to change with time. In that case, the 
scales of the covariance function must be adjusted according to the expected statistical 
evolution of the data field. The choices of scales and variances used in this formulation 
depend on the expected allocation of the field information into signal and noise. In 
this case, the “noise” is actually real information, usually at the smaller scales, that 
is arbitrarily designated as noise. This division depends mostly on what scales are 
present in the data and on what scales we are interested in studying. Additional scales 
can be incorporated explicitly in the covariance function as well, each with a similar 
form to that shown above. The relative weight of each different scale is determined 
by the magnitude of the variance, also known as the amplitude function. Expected 
instrument noise is incorporated in the noise level, E, representing the contribution of 
random noise to the variance at each observation point. The A,, covariance matrix is 
inverted in the mapping process, and a small amount of measurement noise is usually 
needed to keep the process mathematically stable, even if “perfect” data are available. 

An advantage of objective mapping is that it yields a measure of the error, due 
both to the procedure and to the statistical assumptions made about the data. The 
magnitude and structure of this error variance can be used to assess the impact 
of various imposed conditions on the mapping process. Some definitions of terms 
used to describe the accuracy of the mapping are needed at this point. For the 
case of unit signal variance, in which the error maps are non-dimensional, useful 
Mapping is defined as areas with error variance < 0.6 and accurate mapping is defined 
as areas with error variance < 0.3. The signal spatial scale {or scales), Rs, has 
the greatest impact on the estimate field as a whole, both through the need for 


it to be consistent with the actual scales of the data and in its relation to data 
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spacing. The imposed signal scale must be reasonably close to those of the original 
field for the mapping process to work. The area of accurate mapping is considerably 
increased when the signal scale is larger than the data spacing (Bretherton, Davis, and 
Fandry, 1976). Varying the noise scale, Ry, has an impact that seems confined mostly 
to the vicinity of each data point. Choosing a vanishingly small noise scale causes the 
small-scale noise to exert an effect similar to measurement noise over much of the field, 
except exactly at the data locations. At these points it yields a conservative error 
estimate (Bretherton, Davis, and Fandry, 1976). Choosing o% = 0, thereby assuming 
that the “noise” from the small-scale signal is not significant, does not change the 
area of useful mapping, compared to a choice of a small (0% ~ 0.1) noise variance, 
when an infinitesimal noise scale (Ry = 0) is used. However, this vanishing noise 
does increase the area of accurate mapping (Bretherton, Davis, and Fandry, 1976). 
These general considerations were taken into account in the initial selection of the 
covariance function, with the final choice decided by how well the resulting estimate 
field reproduced the original winds. 

To avoid the complication of mapping vector quantities, the U and V components 
were mapped separately as independent scalar properties. Several combinations of 
scales and variances were tested, with an eye towards finding the simplest formulation 
that would accurately reproduce the wind field. The final selection was: 


P=(- 


F(x) = 1.0 exp[—( Sane 








Xx 
500km 


The choice of a 500 km signal scale smooths the original field slightly, depending on 
the amount of variation at yet smaller scales. The choice of o% = 0 reflects the fact 
that we expect the small-scale signal variation to not contribute significantly to the 
total data covariance. After experimenting with different values of the noise variance, 
a simple choice of a constant E = 0.1 was made, so measurement noise is incorporated 
along the diagonal of the covariance matriz, A;;. This value is somewhat arbitrary, 
since it is difficult to precisely determine the measurement, or instrument, noise of 


scatterometer data. Instrument error is a combination of many different SASS sensor 
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data errors, and also includes potential errors introduced during reprocessing of the 
backscatter data. In initial experiments, this choice of a reasonably conservative 
instrument noise variance produced good results, so it was used in all subsequent 
mapping. If the covariances imposed on the estimated field are inconsistent with those 
of the actual data, the mapped field can vary considerably from the original field. More 
severe inconsistency in the imposed conditions can introduce ill-conditioned matrices, 
resulting in a breakdown of the mapping process, either through a failure of the matrix 
inversion, or through the production of worthless estimates. The extent to which the 
estimated field reproduces the original is a good indication of whether the imposed 
scales and variances are appropriate. Ideally, the difference between the estimated 
and original fields at any point should be close to the imposed noise at that point. 
The estimate formed from the simple covariance function and variances given above 
is a good representation of the original data (a more detailed comparison is given 
below), so those choices were employed in all subsequent mapping. 

The chosen mapping procedure is one that includes an estimate of the mean and 
may result in error variances greater than the signal variance (greater than unity in the 
non-dimensional case) in some cases. The mapping process was initially evaluated by 
mapping two revolutions, 744 and 830, that contained a large number of wind vectors, 
and then overlaying vector plots of the original and mapped data and examining the 
structure of the error maps. Overlaying vector plots showed the agreement between 
both the closely colocated vectors and the general circulation patterns of each field. 
In regions of abundant data, the mapping closely reproduceed the original field in 
both respects, with a small amount of smoothing of markedly anomalous vectors, and 
in areas where winds change rapidly over short distances. The latter is particularly 
evident in regions of light and variable winds when no well-defined circulation pattern 
is present. A notable feature of these initial maps is the relatively poor performance 
in data-sparse corners of the estimate grid that require significant extrapolation from 
interior points. These extrapolated border values are noticeably worse than those 
between the swaths, which is an expected result since the extrapolation is one-sided, 


and this emphasizes the need to map only over areas with reasonable amounts of 
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data. Lowest error variance values are found where data is densest, and the double- 
swath pattern of the data field is quite evident in the associated error maps (Figures 7 
and 8). 

Evaluations of this type, made on the initial runs, demonstrated that the chosen 
mapping procedure worked well and that the chosen covariance function was appropri- 
ate. However, the spatial distribution of the data presented an entirely new problem. 
Mapping the diagonal satellite tracks of entire revolutions in the study area produced 
large regions of high error variance, with little useful information, on either side of 
the satellite swath. In addition, the fairly large numbers of wind vectors present in 
each revolution are computationally expensive to map, since the objective mapping 
routine requires inverting a square matrix with a dimension equal to the number of 
data. Thus, the need to discard large amounts of information from the large-area 
Mappings is even less acceptable. The solution chosen here was to map successive 
boxes along the path of each satellite revolution. From one to four boxes were needed 
for each revolution, depending on the length of the track (Figure 9), and data in- 
put for each box included all the vectors within one correlation distance (500 km) of 
each edge. One drawback to this process is that by using smaller grids we lose low 
wave number information, especially in the meridional direction, since the swaths are 
aligned more north-south than east-west. However, this loss should not be significant 
since we expect the smaller, synoptic scale (~ 500 km) to dominate most of the sig- 
nal. Synoptic scale dominance is particularly true in the meridional direction, due 
to the rather strong dependence of weather patterns on latitude and the tendency 
for weather systems to move in a zonal direction. It should be noted that low wave 
number information is also missing in the zonal direction, since the swaths are not 
wide enough to resolve it. The lack of low wave number information in the zonal 
direction cannot be corrected without combining data from several revolutions, and 
in that case, the large temporal and spatial separations involved introduce yet further 
complications. For this study, the loss of low wave number information was accepted 


as a necessary trade-off for the use of a simpler, spatial mapping. 
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As before, the mapped data was checked by comparison to input scatterometer 
data through vector plot overlays and error variance plots, and through comparing 
statistics of the original and estimated fields. This time a more detailed comparison 
was done, using objective maps made for all the revolutions on each of three separate 
days, spread through the data set. The scatterometer revolutions from days 230, 236 
and 248 each contained swaths with large numbers of winds vectors, as well as both 
well-defined circulation patterns and regions of more variable winds. These features 
facilitated the subsequent comparison to the mapped wind field. By overlaying vector 
plots, we can see the agreement between both the closely colocated vectors and the 
general circulation patterns. The circulation patterns, especially those around the 
strong low pressure systems in the north, are closely reproduced, agreeing in the 
inferred location of the low pressure center and associated fronts, and in the direction 
and magnitude of overlapping vectors. Revolution 830, on day 236, is a good example 
of this (Figures 7 and 9). Some of the highest magnitude vectors in the original 
field are not picked up in the estimate field, especially in cases where there is one 
anomalously high vector in a field of generally lower magnitude. Likewise, some of 
the lighter winds are overestimated in the mapped field, particularly in the case when 
they are surrounded by or border on an area of higher magnitude winds. Both of these 
reflect the slight smoothing of the mapping process, which tends to reduce anomalous 
vectors embedded in a more uniform field. More isolated vectors are generally closely 
reproduced in the mapping, since they hold relatively more weight in a data-sparse 
area, and they are then joined by extrapolation to areas with more input information. 
This worked well for vectors that fit reasonably well into the overall pattern and 
did not have extreme spatial separations, which was usually the case for isolated 
points in this data. In any case, the need for excessively large extrapolation was 
minimized, since care was taken to map over regions of denser data. Differences in 
closely colocated vectors from the two fields were most evident in regions of light and 
variable winds, where the winds change rapidly and more randomly over relatively 
short distances. The wind magnitudes tend to be the same, but the directions may 


vary considerably, up to 180 degrees in a few cases. However, since this variation is 
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typical of those areas in the first place, and since the wind magnitudes are small, the 
impact on the final results is minimized. Statistics of the two fields were similar for 
these initial cases, also confirming our choice of covariance matrix. Statistics for each 
field are given for the complete data set in Section 5. Overall, the mapping procedure 
is effective at accurately reproducing the original field; in the general circulation 


patterns, for overlapping vectors, and in a broader, statistical sense. 
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5. Results and discussion 


Prior to the detailed intercomparison of the FNOC and Atlas wind fields, they 
were initially compared to matching synoptic weather charts to obtain a rough subjec- 
tive estimate of their quality. This comparison was used to verify basic flow patterns, 
as well as the location of and circulation around major weather features. It was a 
means to identify gross errors and it also helped to resolve some computer plotting 
questions. To the extent possible, given the coarse nature of the weather maps, both 
wind fields compare well to the flow patterns evident in the weather maps. In particu- 
lar, the series of strong low pressure systems that tracked across the northern section 
of the study area are closely reproduced in the Atlas and FNOC winds. The wind 
field inferred locations of the low pressure centers are sometimes displaced slightly 
(< 200 km) from those of the meteorological analyses, especially for the Atlas winds. 
However, no extreme displacements are noted, and in each case the difference is well 
within the coarse resolution limits of the weather maps. The broader circulation as- 
sociated with these low pressure systems, as reflected in the Atlas and FNOC wind 
patterns, is consistent for each with the isobars and frontal systems of the meteoro- 
logical analyses. This consistency is true both for the spatial patterns observed at 
a particular time, and for the temporal evolution of the patterns across several syn- 
optic analysis periods. The major difference stems from the greater resolution and 
higher variability in the two wind fields, which contain small-scale features that are 
smoothed out of the weather maps. For the Atlas winds, some of the differences can 
also be attributed to the offset in the times of the satellite passes from the 6-hourly 
synoptic analysis times. This difference in time can be up to 3 hours, but is generally 
less than 1 hour for this data. The strongest features during the study period were 
Hurricane Ella and Tropical Storm Flossie, and both showed up clearly and in good 
agreement with the meteorological analyses in the Atlas and FNOC winds. 

The comparison between the mapped scatterometer and FNOC winds was ac- 
complished in several ways, both qualitative and quantitative. Atlas data from a 


particular satellite revolution were matched to colocated FNOC data at the synoptic 
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time closest to the revolution time. Overlaying plots of the two fields provided a good 
qualitative look at their differences. Basically, the fields appear very similar, espe- 
cially across areas with little spatial variation in the wind. Significant differences in 
both magnitude and direction are observed in wind vectors from the two fields in the 
vicinity of strong weather features. Some of these differences can be accounted for by 
the offset in time between the two fields, which may be up to three hours. Also, the 
sampling pattern of the scatterometer did not always provide complete spatial cover- 
age of important features, and filling in with mapped values at the margin of available 
data gives uncertain results, due to the high error variance in the mapping there. The 
difficulty posed by these high variance areas was minimized by mapping only within 
the extent of available data as much as possible. Another explanation of the observed 
differences between the two fields is the higher spatial resolution and therefore greater 
inherent variability of the Atlas winds. Even though there is some smoothing of the 
scatterometer winds, both in the data generation and objective mapping stages, it is 
still on a scale less than that used to generate the FNOC winds. The impact of this 
greater variability is expected to be much more near areas of intense, smaller-scale 
circulation, and the observed differences in the two wind fields reflect this. 
Quantitative comparison of the Atlas and FNOC wind vectors is accomplished 
in several ways: 1) statistical descriptions of each field, 2) differences in the mean 
values within latitude-longitude boxes, 3) scatter plots of wind speed and direction for 
each field, and 4) statistical descriptions of the difference field (Atlas-FNOC). Initial 
evaluations include the ungridded as well as the mapped Atlas winds, but subsequent 
calculations require colocated vectors and are only done for the mapped Atlas and 
FNOC winds. Most calculations were done on the east and north components of the 
wind, but some vector statistics are also shown. The various quantities are computed 
for the instantaneous vector field “snapshots” and also for one-day and three-day 
averaging periods. In addition, the statistics are also subsequently grouped by latitude 
to examine its impact, since weather patterns, and therefore wind, tend to group in 
zonal bands. Most of the comparisons done in this study concentrate on the wind 


velocity, since that is the data available, rather than on the wind stress, though it is 
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recognized that stress is a better representation of the driving force provided by the 
wind. However, the conversion from velocity to stress is not trivial, especially over 
larger areas and longer periods of time, and introduces additional uncertainties due to 
factors such the formulation of the drag coefficient. Also, the stress values are quite 
small and this can cause problems in the difference calculations. For these reasons, 
the bulk of the statistics are calculated for the wind velocity, either in component 
or vector form. In addition to these velocity calculations, the wind stress curl is 
calculated for the Atlas and FNOC wind fields, and examined for significant and 
consistent differences betweem them. Since this parameter is so crucial to ocean 
circulation forcing, we found it useful to find and compare the values from the two 


fields, recognizing the limitations on stress calculations discussed above. 
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Field statistics 


Three separate wind fields are examined in this study, the ungridded Atlas SASS 
winds, the objectively mapped Atlas winds (derived from the first), and the FNOC 
combined pressure-derived and buoy and ship report winds. Descriptions of the indi- 
vidual wind fields are provided in various forms, qualitative and quantitative. Repre- 
sentative vector plots, mean and standard deviation for data within individual areas 
(matching the objective mapping grid), and the mean and standard deviation of all 
data on a particular day are given for all three fields. In addition, for the mapped 
Atlas and FNOC fields, mean speed and vector statistics, daily averages grouped by 


latitude, and three-day averages are shown. 
Ungridded Atlas winds 


For the ungridded Atlas winds, Table 1 shows the U and V component statistics 
for individual areas in each revolution, including mean, variance, weighted variance, 
standard deviation, and number of data. In this and following tables, the variance 
is simply the sample variance and the weighted variance accounts for the effect of 


varying wind speeds on variance values and is normalized by the rms value: 


The data used includes all those vectors within the corresponding mapping grid as 
well as within one correlation length (that used in the mapping) of any grid edge. This 
is the same data that was used as input to the objective mapping program. Figure 10 
shows the distribution of the component means and standard deviations graphically. 
Areas of high wind speed, such as 715 U1 and 735 U1 (where U1 indicates a particular 
box along the satellite revolution), or high variance relative to wind speed, such as 
715 V1 and 779 U2, can be associated with meteorological features present near those 
areas, as shown by synoptic weather maps. Those areas with high mean wind speed 
and significant variance, especially 836 U2 and V2 and 1008 U2 and V2, straddle 


circulation about well-developed low pressure systems, reflecting the strength of the 
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system and the sharp variation in wind direction across fronts and around the low. 
Areas that include a single front, like 715 1, show a large mean and lower variance in 
one direction (U in this case), and a low mean with much larger variance in the other 
direction. Across this particular front, while the winds remain predominantly easterly, 
the north-south direction varies up to 180°, driving down the mean and increasing 
the variance in the V component. For the entire data set, the mean U values tend 
to be greater than the corresponding V values on average, roughly two-thirds of the 
time, but the weighted variances are split nearly evenly in relative size for the two 
components. 

Statistics for the U and V components of the ungridded Atlas winds, averaged 
over one-day periods rather than just instantaneous “snapshots”, are listed in Table 2. 
The mean magnitudes are evenly split, with neither U or V relatively larger on average, 
and the weighted variances are fairly even, with 02 > a? 17 out of 29 times (59%). The 
daily means are lower than those for the individual areas in Table 1, and the variances 
increase, with the standard deviation now greater than the mean in each case. The 
lower mean values and higher variances are consistent with what we expect, since over 
the course of one day the samples are taken from widely spaced areas, usually from 2 


or 3 revolutions, and the winds can vary greatly within those space and time scales. 
Mapped Atlas winds 


The distribution of U and V components, and wind speed and direction from all 
the mapped Atlas vectors in this data set are shown in Figure 11. The U component is 
evenly divided between positive (westerly wind) and negative, while the V component 
is more frequently positive (southerly wind). The distribution of magnitudes coarsely 
reflects the Rayleigh distribution expected of the original Atlas data, providing more 
confirmation of the objective mapping process. The wind direction bins indicate the 
compass direction the wind is blowing toward, for example, 0° is a southerly (or 
northward) wind. This distribution also shows the slightly more frequent occurrence 


of southerly winds. Various statistics for the U and V components in individual 
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boxes (within revolutions) for the objectively mapped Atlas SASS winds are listed in 
Table 3 and shown graphically in Figure 12. It should be noted that these cannot be 
directly compared to the corresponding values for the ungridded winds, since those 
were calculated over slightly larger areas and included more data. Again, mean U 
magnitudes tend to be larger, with | U |>| V | 68 out of 117 times (59%). The relative 
values of the weighted component variances are generally close, and o7 > o2 51 out 
of 117 times (44%). Compared to the values for the ungridded winds, the weighted 
variances of the mapped winds tend to be slightly lower, with o2, > Gaus 187 out 
of 234 times (80%), including 101 of 117 (86%) for U and 86 of 117 (74%) for V. 
This trend is expected since the objective mapping process smooths the data field 
to some extent. Table 4 presents the statistics for the U and V components of the 
mapped SASS winds, averaged over one day periods. The means are reduced relative 
to the individual “snapshot” values, and the variances increase to the point where 
the standard deviation is greater than the mean in most, though not all, cases. The 
component magnitudes are evenly split in relative size, while the weighted variances 
show some difference between components, with o2 > a2 19 out of 29 times (66%). 
The greater variance of the zonal wind, on average, reflects the slightly greater range of 
speeds associated with it, as compared to the meridional wind. In North Atlantic mid- 
latitudes, one expects the average zonal wind to be larger than the average meridional 
wind (Leetmaa and Bunker, 1978). Component statistics, again averaged over one day 
periods and now also grouped into 10° latitude bins are shown in Table 5. Tables 6 
and 7 list one-day and three-day average vector statistics, including the mean speed, 
the speed and vector variance and standard deviation, and the percentage of the 
vectors that fall in each of four direction bins. The standard deviation of the wind 


vector is defined as: 
ly = =) 2)\ + 
Ovec = (w—7 du = Wey (gee |) 
while the standard deviation of the mean wind speed is defined as: 
eo = ed 
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where: | U |= (u? + v?)3 


The mean speed and vector deviations are quite different, with o,,q smaller and 
much more uniform across a given area than oy... In this case the directions are 
given in the oceanographic sense of going toward a particular direction. The bins 
each include 90° of arc, centered on the given direction. For example, the north bin 


includes 315° < @ < 45°, where @ is the direction that the wind is blowing toward. 


FNOC winds 


The distributions of the U and V components, and wind speed and direction from 
all the FNOC wind vectors are shown in Figure 13. As with the Atlas winds, the U 
component distribution is divided evenly between positive and negative values and the 
V component is more frequently positive, indicating southerly winds. The magnitude 
distribution differs from that of the Atlas winds, with more winds at the lower wind 
speeds. However, there are several FNOC values in the highest speed range (25-30 
m/s), while there are no Atlas values of that magnitude. The direction distribution 
shows the tendency for the FNOC winds to be southerly, to an even greater extent 
than the Atlas winds. Statistics for the U and V components in individual areas for 
the FNOC wind field are listed in Table 8, and shown graphically in Figure 14. The 
chosen times are those closest to the time of each satellite pass, and given the 6-hour 
spacing of the FNOC analyses, may be up to 3 hours different from the time of the 
satellite pass. Most times were much closer than this upper limit, and in all cases, 
the best FNOC match was at either the 00Z or 12Z analysis times. The mean U 
magnitude is larger than the mean V magnitude slightly more times, 88 of 144 (61%), 
with one case where both are the same, and a slightly greater proportion of the V 
component winds have higher variances than the corresponding U component, 83 of 
145 (57%). The trend observed in the two sets of Atlas winds, of higher means in 
one component relative to the other being associated with lower relative variances 


in that same component, is also seen in the FNOC data. Again, we expect winds 
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of greater variability to have their means reduced, when cases of both positive and 
negative component wind speeds are summed. 

Table 9 presents the one-day average statistics for the U and V components of 
the FNOC winds. As was true for the scatterometer winds, the daily means are lower 
than the earlier means for individual areas, and the variances are relatively larger, 
with standard deviations greater than the mean values in most cases. The mean 
U value tends to be larger than the corresponding V mean value more of the time, 
with |U| > |V|, 17 of 29 times (59%). The weighted variance of the U component is 
also greater than that of the V component more frequently, with a? >a? 19 of 29 
times (66%). The greater variance in the zonal wind can again be associated with 
the wider speed range of the zonal (-2.5 to 3.7 m/s) compared to the meridional (-0.8 
to 3.8 m/s) winds. The mean daily meridional winds are mostly positive, indicating 
southerly winds, which is the expected average flow in this area at this time of year. 
This behavior parallels what is observed in the scatterometer winds as well. Table 10 
lists statistics for the component winds, averaged over one day periods and grouped 
into 10° latitude bins. Tables 11 and 12 list one-day and three-day average vector 
statistics, including the percentage of the vectors within each averaging period that 
fall in each of four direction bins. Again, the directions indicate wind blowing toward 


a particular direction. 
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Boxed mean differences 


Prior to regridding the scatterometer winds, a preliminary evaluation of the dif- 
ference between the scatterometer and FNOC winds was made by calculating the 
difference between the mean U and V components over 10° latitude by 10° longitude 
boxes. This was done for a limited subsection of the data, three separate days, in 
order to get a rough idea of the magnitude of the difference and possibly any obvious 
bias or structure. To increase the amount of data used in calculating the mean values 
and to provide adequate spatial coverage, daily averages were computed. This was 
done simply by using all the FNOC data from the 00Z and 12Z, as well as 00Z of 
the following day, synoptic analysis times, and the scatterometer data that fell within 
three hours of those times. Data from three separate days, 18 and 24 August and 
5 September (Julian days 230, 236, and 248 respectively), were investigated. For 
day 230, data from SEASAT revolutions 744, 750, 751, and 758 were used; for day 
236, from revolutions 830, 836, 837, and 844; and for day 248, from revolutions 1002, 
1003, 1008, 1009, 1016. Tables 13, 14, and 15 list the AU and AV mean and stan- 
dard deviation for individual latitude-longitude boxes and the entire area. The AU 
(U, — U;) and AV (V, — V¢) for each box are also displayed in a grid. Some boxes 
contained no scatterometer values for a particular day and are marked N/A. For days 
230, 236, and 248, the respective AU values are -1.29, 1.12, and -1.03 m/s, while the 
AV values are 0.57, 1.34, and 2.14 m/s. The average values of the boxed differences 
are somewhat lower, with —0.63 + 2.04, 0.36 + 1.19, and —0.91 + 1.93 m/s for mean 
AU, and 0.05 + 2.09, 0.58 + 2.13, and 0.52 + 2.44 m/s for mean AV, with all mean 
differences less than one standard deviation. 

The difference between mean values for the entire region is thus larger than 
the difference expected within the smaller area of any latitude-longitude box. The 
increased difference is due in part to the greater spatial variability of the wind field 
when averaged across the larger area. The average component variances computed 
for the entire area, for both SASS and FNOC winds, are close in magnitude and 


day-to-day differences, with the only significant difference on day 248, between the 
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U,, variance of 46.23 m?/s? and the U; variance of 29.98 m?/s*. The boxed mean 
differences exhibit no obvious latitudinal, meridional, or temporal biases or trends, at 
least not evident over these three days. The boxed variances of the Atlas winds tend 
to be higher than the corresponding FNOC wind variances, for both U (23 out of 30) 
and V (22 out of 30) components. The relatively greater variances of the Atlas winds 
is an expected result, since the higher resolution scatterometer winds should be more 
variable than those from the smoothed FNOC surface analysis. Those cases where 
the FNOC variances are higher generally correspond to coverage of strong circulation 
features, picked up completely by the FNOC analysis on its regular grid, but only 
partially covered by the satellite swath. This potential, but unresolvable, difference 
in spatial coverage within a given latitude-longitude box is one reason for regridding 
the satellite data to locations corresponding to those of the FNOC data. From this 
very limited set of data the most significant fact to emerge is the negative result 
that no readily apparent bias exists between the two fields. The limited methods 
of intercomparison (ie: boxed mean values) and associated poor resolution of the 
difference field dictate an alternative process, in which scatterometer and FNOC winds 


are colocated on a regular grid to facilitate statistical comparisons. 
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Scatter plots of wind speed and direction 


As a check on the mapped winds, the average daily standard deviations of the 
U and V components of the original and mapped scatterometer winds were plotted 
versus one another (Figure 15). Values fall along the 45° line without much scatter, 
providing confirmation of the chosen mapping procedure. Scatter plots of Atlas vs 
FNOC wind speeds and directions over one-day periods are given in Figures 16 to 
44. The speed plots simply show the degree of agreement about a 45° line, while the 
direction plots have a slightly more difficult interpretation, since the periodic nature 
of compass directions means that values in the upper left and lower right corners 
also indicate closer agreement. The most obvious feature of nearly all the plots, 
especially those for direction, is the large degree of scatter present. The scatter tends 
to increase at higher wind speeds, in some cases about the 45° line, and in others 
tailing off to one side. With the exception of days 229, 230, 231, 244, 247, and 248, 
the daily averaged Atlas wind speeds are greater than the FNOC wind speeds. No 
consistent trends are evident in the pattern of agreement as a function of wind speed. 
On some days (see days 225 and 235) the scatter seems evenly split between Atlas 
and FNOC, with or without offset from the 45° line. On others (see days 232, 233, 
and 234) the pattern trails to the right at higher wind speeds, indicating that the 
higher Atlas winds were not picked up by the FNOC analysis. Day 248 shows an 
especially distinctive curve to the right for the highest Atlas wind speeds. During this 
period Hurricane Ella was moving along the east coast of North America, weakening 
and moving northeast over the open ocean by 12Z on day 248. The Atlas winds at 
this time, from revolution 1008 (Figure 45), are considerably higher than those from 
the corresponding FNOC analysis (Figure 46). The offset in time should not be a 
significant factor, since this satellite pass was only ~ 40 minutes prior to the analysis 
time. The difference between the two fields can be interpreted as either a too rapid 
weakening in the meteorological analysis, or an overestimation by the scatterometer, 
perhaps due to the extreme surface weather associated with the hurricane. Without 


additional information the exact cause cannot be determined with certainty, but it 


49 





is a good example of possible differences in the two fields in response to a strong 
circulation feature. 

There are also several days (see days 230, 246, and 247) in which high FNOC 
winds are not reflected by Atlas values, tilting the scattered points to the left. Day 244 
shows this trend of greater FNOC speeds, at higher wind speed ranges, in addition 
to a very large scatter at the higher wind speeds. For day 230, revolution 744 in 
particular, the tilt seems to be due to the mapping, more than to a difference in the 
actual winds. The most northern box was extrapolated too far at its northern edge, 
and since the wind field magnitudes increase to the north, the extrapolated estimates 
are significantly lighter than the corresponding FNOC winds. The objective map of 
revolution 744 was one of the first made, and demonstrates the artificial differences 
that can be imposed by even slightly faulty mapping. More care was taken with 
subsequent maps to keep the boundaries of the estimate grids within the extent of 
the input data. For days 244 and 247, the likely cause of the higher FNOC winds 
is again related to the passage of Hurricane Ella. On day 244 the hurricane was 
intense, and small in area, and the rapid change in winds over short distances could 
easily account for the large scatter seen in the Atlas vs. FNOC scatter plot at the 
higher wind speeds. Strongest winds were present during this period in the fully 
developed hurricane, and limited comparisons to high-quality surface wind data have 
shown some tendency for the SASS to underestimate the highest (> 20 m/s) wind 
speeds (Jones et al., 1982). This may be the case in this situation, especially since, as 
shown above, when the winds abated by day 248 the Atlas values became relatively 
larger than those of the FNOC analysis. For day 246 (with no satellite passes over 
Hurricane Ella), the probable cause of the relatively high FNOC winds is that the 
strongest winds in a developing tropical storm at the southern edge of the grid area 
were not picked up in the available Atlas data. The consistent differences exhibited by 
these specific cases were not reflected through the entire data set. Despite the large 
amount of scatter about any estimated regression line on a particular day and the 
different slopes of those lines for different days, the overall trend is towards agreement 


between Atlas and FNOC wind speeds across the entire data set. 
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Agreement in wind directions also follows the same general trend, with even more 
variation and no easily discernible pattern. Differences of up to 180° are found, es- 
pecially on days 228, 229, 230, 234, 235, and 251. Two factors can strongly influence 
directional differences, both due to slight mismatches in the sampling times and areas 
of the two fields. The first is the inherent variation under light and variable wind 
conditions (under high pressure cells for example), in which small differences in sam- 
pling times and locations, and integrating intervals could easily show up as observed 
differences. The second factor involves offsets in the locations of well-developed circu- 
lation patterns for the two fields, especially cyclones, where the actual winds sharply 
change direction over relatively short distances. The potential for differences from 
offset patterns is particularly large for intense features such as hurricanes, where the 
change of wind direction is extremely rapid and the wind speeds are very large. The 
differences resulting from offset patterns are “real”, since the winds seen by each field 
are not quite the same. For most of this data set, the meteorological features tend to 
match up well, so the first factor has a more significant role in the average directional 
differences. The exception is again related to the passage of Hurricane Ella, and offsets 
in its location in the two fields lead to extreme directional differences. The difference 
in circulation patterns is evident in the Atlas winds from revolution 1008 (Figure 45), 
where the inferred location of the low pressure center is to the northeast of that of the 
day 248 12Z FNOC analysis (Figure 46). The broader circulation patterns associated 
with less intense cyclones are generally much more similar between the two fields and 
exhibit less extreme differences than the hurricane case. Overall, the comparisons of 
the individual wind vectors shown in these scatter plots show that there are significant 
differences in the fields, particularly in the generally larger magnitude of the Atlas 
winds. The wide scatter prevents resolution of other trends, however, at least of any 
pattern that remains consistent throughout most of the data. 

In addition to the daily plots of wind speed and direction comparisons, various 
averaged quantities were plotted in the same fashion. These averaged values are shown 
in Figures 47 to 58, and include various one-day and three-day period plots of mean 


speed and speed and vector standard deviation, as well as the means and standard 
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deviations for the U and V components, including a separation into latitude bins. The 
form of the mean speed and vector standard deviations are the same as given for the 
field statistics. From the daily mean speed plot (Figure 47) it is obvious that the Atlas 
winds are higher on average. The spatial variability, as expressed by the standard 
deviation of the mean speed, seems to be generally the same, given the small sample 
size and fair amount of scatter (Figure 48). The vector standard deviation shows that 
the Atlas wind directions are slightly more variable on average than the FNOC values 
(Figure 48). The three-day average values show a similar pattern, with all average 
Atlas speeds greater than FNOC speeds, a nearly even split in the standard deviations 
of the mean speed, and all but one of the Atlas vector standard deviations larger than 
the corresponding FNOC value (Figures 49 and 50). Scatter plots of the daily mean 
values of the U and V wind components (Figure 51) show a trend of agreement along 
the 45° line for U and a wider dispersion with more highly negative Atlas than FNOC 
values for V. For both components, and most markedly for eastward velocity, most 
of the average FNOC values are larger than the corresponding Atlas velocities. This 
can be partially attributed to the greater variability of the Atlas winds, however, 
which tends to reduce the mean U and V values as negative and positive velocities 
cancel each other. Eastward average velocities are split rather evenly between positive 
and negative, while northward velocities are mostly positive. This corresponds with 
expected average flow patterns over this area at this time of year. 

Plots of the standard deviations of the U and V components (Figure 52) show 
general agreement between Atlas and FNOC variability, with greater dispersion for 
the northward values. For both components, especially the northward, most Atlas 
values are greater than the FNOC values. The relatively greater variability of the 
Atlas winds is expected, since the FNOC fields are smoothed over a larger area than 
the Atlas fields. The plots of U and V component means and standard deviations, 
grouped into latitude bins (Figures 53 to 58) exhibit many of the same patterns as 
the combined daily values. The Atlas standard deviations are greater and the mean 
FNOC values are slightly higher. The eastward velocities are negative (southerly) 


in the 20-30°N bin, increase and become positive in the 30-40°N bin, and decrease 
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slightly in magnitude and scatter about zero in the 40-50°N bin, averaging positive in 
each latitude bin. The standard deviations for both components are smallest in the 


20-30°N bin and of relatively similar magnitude in the two northern bins. 
Difference field statistics 


Various statistics on the U and V components of the Atlas-F NOC vector differ- 
ence field for individual areas within each revolution are listed in Table 16, including 
the mean, variance, weighted variance, standard deviation, mean squared difference, 
and weighted mean squared difference. The distributions of the AU and AV com- 
ponent means, and the speed and direction of the difference vectors (formed from 
individual AU and AV components) are shown graphically in Figure 59. Most values 
of AU, AV, and IAU| fall within a +5 m/s range, with a much smaller, but signifi- 
cant number in the -10 to -5 m/s range for AV and 5 to 10 m/s range for |AU|. The 
difference vector direction, 9,4,, is evenly distributed among the various direction 
bins, with slightly more in the 180° — 300° range (southeasterly in the meteorological 
sense). So the difference vector does not appear to have any preferential direction over 
this data set. The revolution and box number of the Atlas wind field and the synoptic 
analysis time of the matching FNOC wind field are listed in each case. The mean 
squared difference (MSD) is simply the expected value of the square of the difference 
between either the U or V Atlas and FNOC components: 


N 


1 2 
MSD = ~ > (Uai — Uri) 


1 
The weighted mean squared difference (WMSD) is weighted by the square of the Atlas 


field: 
Ye (Uai — Uri)? 
yn 


As with the weighted variance, the normalized mean squared difference is more useful 


WMSD = 


for intercomparisons of values from different areas, since it removes the influence of 
wind strength on the MSD magnitude. Areas of stronger winds will tend to have larger 


mean squared differences than those of lighter winds, simply because the individual 
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magnitudes are larger, but a more accurate measure is to determine which of two cases 
is proportionally larger, relative to the respective wind strengths, rather than in the 
absolute sense. The relative sizes of the means and variances of the two components 
are evenly split, with |U| > |V| 71 out of 144 times (49%) and o2 > o2 70 out of 141 
times (50%), with three cases of equal variance. 

One-day averaged statistics for the U and V components of the vector difference 
field are listed in Table 17. The mean daily U and V magnitudes are quite small 
and less than one standard deviation in nearly all cases, with the only exceptions 
coming in the V component means on days 237 and 243. As expected, the one- 
day average means are lower than the “snapshot” values from individual areas, and 
the corresponding variances are higher than those for individual areas. The large 
variances of the averaged difference field components, with most weighted variances 
greater than 0.9, combined with the low mean values, indicates a lack of obvious 
bias between the two fields over one day periods. The values of the U component 
range from -1.73 to 1.15 m/s and those of the V component from -3.13 to 1.36 m/s, 
with the largest magnitudes occuring during days 242 to 245, particularly in the V 
component. Though the V component has a slightly greater range than U, neither 
component is consistently larger than the other, with |U| > |V| 13 out of 29 times 
(45%). The relative sizes of the weighted variances are also split evenly between the 
two components, with o? > o? 15 out of 29 times (52%). The mean U component 
is nearly evenly split between positive (westerly winds) and negative (easterly winds) 
directions, with U positive 13 out of 29 times (45%), while the mean V component 
shows more bias in sign, with V positive (southerly winds) only 8 of 29 times (28%). 
So, though there does not seem to be an east-west directional bias between the Atlas 
and FNOC winds, the Atlas winds tend to be more southerly (negative) than the 
FNOC winds. Table 18 presents U and V component statistics, grouped into 10° 
latitude bins and averaged over one-day periods. 

In addition to various statistical quantities calculated for the U and V components 
of the difference field, mean speed and vector statistics were also computed. To 


calculate these values, the component difference field (AU and AV) was treated as a 


o4 





vector, with the speed, AU, calculated as (AU? + AV} a. and the speed and vector 
variances calculated as shown in the preceding section. This method of calculating 
the mean speed generally yields values considerably larger than the difference of the 
magnitude of the Atlas and FNOC vectors, especially in cases where the Atlas and 
FNOC components are of opposite sign. So these mean speed vaules reflect differences 
is direction as well as magnitude. Table 19 contains these values for individual areas 
and instantaneous times, Table 20 has the same for the entire area averaged over 
one-day periods, and Table 21 contains values grouped into 10° latitude bins and 
averaged over one day periods. The largest mean speed values for individual areas, 
such as 952 2, 988 1, 1008 2, and 1074 2, are all associated with circulation around 
strong low pressure systems, with Hurricane Ella in the first three cases. The rapid 
change in direction of these winds, across relatively short distances, results in large 
differences if the location of the low pressure center is even slightly offset between the 
Atlas and FNOC fields. A good example of this is Atlas revolution 1008 (Figure 45), 
mentioned in the previous section, where the inferred location of the low is displaced 
to the northeast, relative to that of the FNOC day 248 12Z analysis (Figure 46). The 
standard deviation of the mean speed is always less than the mean value, with the 
largest variance in the areas of highest mean speed, such as those discussed above. The 
very large vector variance values for all these same areas also establish the significant 
impact of the directional differences. The one-day averages do not exhibit the extreme 
values of the individual areas, with the average mean speed remaining the same, and 
the speed and vector variances increasing slightly on average. The standard deviation 
of the mean speed is less thant the mean speed in all but one case, day 248, which 
is associated with the passage of Hurricane Ella. An interesting trend in the daily 
averages is that both the speed and vector variance increase towards the end of the 
study period (late August and early September). The progressive increase is most 
likely associated with the stronger weather patterns that begin to track further south 
at that time (reaching the northern edge of the study area), as well as with the 
influence of Hurricane Ella. The daily averages, grouped by latitude, exhibit similar 


trends, and the standard deviation of the mean speed is greater than the mean value 
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in every case. The most noticeable trend in these is the distinct increase of the mean 
speed, and speed and vector variances, from south to north. The largest values fall 
within the 40° — 50° N bin, and this supports the association of greater Atlas-FNOC 


differences with the stronger weather patterns that are generally found there. 
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Wind stress cur] 


A primary use for synoptic scale wind fields is as a driving force for ocean circu- 
lation models, since the wind stress is the single most important factor in determining 
that circulation. The input of wind energy into the ocean through the surface stress 
is reflected in the directly wind-driven Ekman flow, in the geostrophic flow due to 
divergence of the Ekman flow, and in the indirectly wind-driven motions due to ther- 
mohaline forcing. Considering just the frictional Ekman layer at the surface, vertical 
integration of the Ekman equations, in combination with mass conservation, yields a 
simple relation between the vertical component of the curl of the wind stress and the 


vertical velocity at the base of the Ekman layer: 
OM keurl? 


where: 
f = Coriolis parameter (1/s) 
p = water density (kg/m°) 
We = vertical velocity at base of the Ekman layer (m/s) 


kcurl7 = vertical component of the wind stress curl (N/m) 


The vertical velocity at the base of the Ekman layer, driven by the wind stress curl, 
then acts as an upper boundary condition on the interior circulation. In the interior 
of the ocean, beneath the Ekman layer, we expect the linear vorticity balance to hold: 


DN is 


where: 
( = variation of the Coriolis parameter with latitude (1/ms) 
v = meridional velocity (m/s) 


Ww : : ° : 
—— = vertical gradient of vertical velocity 


Oz 


In the general case, the wind stress curl acts as a driving force on this interior flow, 


through the dynamics of the Ekman layer, but a special case illustrates a more direct 
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connection between the circulation and wind stress curl. The Sverdrup relation is a 
simple form of the relation between the curl of the wind stress and the mass transport, 
found by vertically integrating the linear vorticity equation, neglecting lateral friction, 
and assuming that vertical velocity vanishes at the bottom. The resulting equation 
is: 


~ 


BMy = keurl7 


where: 


M, = vertically integrated meridional transport 


As discussed by Stommel (1965), the Sverdrup relation actually expresses a balance 
between the divergence of the Ekman transport and the divergence of the geostrophic 
transport. A further discussion of this dynamical balance and the relationship be- 
tween the wind stress curl, Ekman transport, and geostrophic transport is given by 
Stommel (1965). It must be noted that, in addition to the frictional wind stress-driven 
motions described by these relations, motions forced by thermohaline processes are 
also significant. These are not discussed further, except to note that the surface wind 
field plays an important indirect role in this forcing as well, through its impact on 
the surface fluxes and on mixing in the water column. Further discussion of some 
basic aspects of thermohaline processes can be found in Pond and Pickard (1983) and 
Veronis (1981). 

These relationships are then used by modelers to impose the direct and indirect 
forcing of the wind stress as a boundary condition on the ocean circulation. Pond 
and Pickard (1983) provide a brief review of various aspects of this process, including 
the typical wind fields used, the impact of wind field resolution, and the wind stress 
forcing used in several models. The usual choice of wind stress information is a cli- 
matological average wind stress field, several of which are available (Bunker, 1976; 
Han and Lee, 1983; Hellerman and Rosenstein, 1983). The shortest time scale of 
these fields is usually a monthly average, which is useful for computing mean circu- 
lation patterns, but which has limited application at shorter time scales. Large-scale 


ocean currents and ocean properties such as the sea surface temperature and vertical 
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temperature distribution display pronounced spatial and temporal variability. The 
variations of ocean motions and properties are directly related to the variability of 
the surface wind field, most of which is lost in the averaged climatological fields. The 
coarse spatial and temporal resolution associated with climatological fields prevents 
their use in examining the detailed structure and evolution of the ocean circulation in 
response to wind forcing. In addition to this resolution problem, even the accuracy 
of climatology is questionable, since it does not necessarily reflect local conditions 
at any particular time. The results of forcing from an averaged wind field may be 
considerably different from the overall effect of a temporally evolving wind field. The 
synoptic fields discussed in this paper are more appropriate than the climatological 
mean fields for forcing on a daily basis, as is needed to update numerical forecast 
models for example. Since knowledge of the wind stress curl is so crucial, it is impor- 
tant to examine the wind stress curl values for the Atlas and FNOC wind fields, in 
addition to the wind velocities. 

For each field, the first step is to calculate wind stress from the wind velocity, 


using the wind stress and drag coefficient formulations of Large and Pond (1981). 


i pCp|U|U 
where: 
ee pCp|Ulu 
Ty = pCp|Ulv 


U = vector wind velocity 


Cp = drag coefficient based on total wind magnitude, |U| 


The wind magnitude, |U|, acts as a scale factor for each of the stress components, 
reflecting the fact that each velocity component will have an effect on the surface 
roughness, and that the combined effect will appear in both stress components. For 


the same reason, the drag coefficient, Cp, as a measure of the surface roughness, is 
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based on the magnitude of the wind, rather than the individual U and V components. 


Once the stress is found at each point, the curl is calculated as: 


OT, OTR 


where the partial derivatives are estimated by a simple finite difference. Using a 
finite difference method may underestimate the wind stress curl, especially if there 
are large spatial separations in the data. The underestimation is‘more of a concern in 
curl calculations made from wind stress averaged over longer periods, which tend to 
be more widely spaced than the 1° latitude by 1° longitude grid used here. Though 
the amount of underestimation due to finite differencing is reduced for the smaller 
grid spacing, it will still be present to some degree, yielding slightly conservative wind 
stress curl values. 

With this scheme, the wind stress curl was computed at the interior points of each 
individual grid area for the Atlas and FNOC winds. These were then averaged for 
each grid and also over one day periods. Values for the wind stress curl, averaged over 
individual areas for the Atlas winds, are listed in Table 22. Included are the mean, 
variance, weighted variance, standard deviation, and the number of data points. These 
statistics are calculated in the same manner used previously for the wind velocity 
components. The average curlr4 magnitude of 1.2 x 107-’N/m° is in agreement with 
the O(10~’) magnitudes generally expected (Pond and Pickard, 1983). The wind 
stress curl is positive in 96 of 145 times (66%), which is an expected result, since 
estimates made from climatological winds tend to show positive wind stress curl north 
of approximately 30°N in the North Atlantic (Leetmaa and Bunker, 1978). The mean 
values are generally close to the respective standard deviations, with the standard 
deviations being larger a slightly greater percentage of the time, og > |curlrg| 85 of 
145 times (59%), with both the same three times. 

Wind stress curl averages for individual areas of the FNOC winds are shown in 
Table 23, in the same format as for the Atlas values. The average curlrp magnitude of 
1.0 x 10-’N/m° is again of the expected order of magnitude, and it is positive 103 of 


145 times (71%). As with the Atlas values, the mean curlrp values are generally close 
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to their respective standard deviations, with op > |curl7p| 74 of 145 times (51%), and 
equal three times. For both the Atlas and FNOC winds, the relatively large standard 
deviations indicate that there is considerable variation of the wind stress curl, even 
within the gridded areas. Comparing the relative magnitudes of values from the two 
fields in each area, |curl7a| > |curl7p| 91 of 145 times (63%), with both the same 
three times. Both in an average sense (from all areas) and for individual grids, the 
Atlas winds tend to have aslightly greater wind stress curl magnitude than the FNOC 
winds. In general, the respective area averages from each field agree in sign, with both 
having the same sign 122 of 145 times (84%). The average difference (in magnitudes) 
between mean curl values of the two fields is 8.1 x 107-8N/m?®, which is of the same 
order as the means, but is also less than the standard deviations, so the difference is 
often swamped by the variation of the two fields. 

The Atlas and FNOC wind stress curl fields were also compared by examining 
contour plots of the wind stress curl values within individual boxes. For both the 
Atlas and FNOC fields, results from the contour plots are consistent with the area 
statistics, also reflecting the increased magnitudes and increased spatial variability of 
the wind stress curl in regions of more intense circulation. More quiescent areas are 
associated with broad patterns and smaller curl magnitudes, such as revolution 664 2, 
while areas of stronger weather are associated with more variable patterns and larger 
curlr magnitudes, such as 664 3 (Figure 60) or the extreme case of 715 1 (Figure 61). 
There is considerable variation in the amount of agreement between the plots for each 
field, with no obvious, consistent pattern to the differences. There is generally better 
agreement in regions of broad flow, with little spatial variability (revolution 1059 2 
- Figure 62), partially because the Atlas and FNOC wind fields usually correspond 
well in those cases. Frequently, the basic structure of the two curl7 fields is similar, 
especially in the latitudinal variation, but colocated individual values can be quite 
different, either because the range of values is different or there is a bias between the 
two fields. 

Revolution 830 1 (Figure 63) illustrates both the similarities and differences that 


can be present in the maiching contour plots. In this case, the same general trends are 
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reflected in each field, including similar magnitudes across most of the area, the same 
tendency for curlr values to group in meridional bands, and similar zonal variation of 
the sign of the values. Large negative values are present in the east, becoming more 
positive to the west, and then again more negative yet further west. However, distinct 
differences are also present along the western edge, in the more negative Atlas values 
and the greater meridional variation of the FNOC winds. Also, the meridional bands 
do not have quite the same orientation for each field, and the values of the two fields 
are slightly different, such as along 52°W, where the postive Atlas values reach much 
farther south than the positive FNOC values. A similar case is revolution 1008 2 
(Figure 64), where the two fields have similar patterns, but the FNOC values seem 
displaced to the southeast relative to the Atlas values, and the Atlas values are more 
extreme, especially near 45°N 48°W and in the northeast corner. Despite the overall 
correspondence between the Atlas and FNOC curlz fields seen in many plots, there 
are also many areas (both with light and strong winds) where neither the patterns 
nor values of the two fields are in agreement, such as revolution 1008 1 (Figure 65). A 
consistent difference between the two fields is not apparent in the curly contour plots, 
so a general conclusion cannot be drawn. 

The point values of wind stress curl were also averaged over one-day periods, and 
the results are listed in Table 24 for the Atlas winds and in Table 25 for the FNOC 
winds. The same statistics are computed as in the previous case. As for the area 
averages, the daily averages tend to be positive, with curlr, and curlrp each positive 
19 of 29 times (66%). For both the Atlas and FNOC data, the standard deviation 
is greater than the mean for every day. Averaging the magnitudes of all the daily 
values also shows that the standard deviation is greater than the mean, with loal = 
2.9x 10-7N/m? compared to |curlra| = 6.7 x 10-°N/m’, and |op| = 3.3 x 1077N/m? 
compared to |curl7p| = 5.3 x 10-®N/m*. This large variation is most likely due to the 
wide spatial separation that occurs over the course of one day as data from different 
revolutions are combined. The relative magnitudes of the daily average wind stress 
curl reflect the same trend present in the individual areas, with |curlr,| > |curl7p| 19 


of 29 times (66%). The average |curlr,| value is also slightly larger than the average 
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lcurl7p| value so, as for the individual areas, the Atlas wind stress curl tends to be 
slightly greater. The daily averages from each agree in sign most of the time, with both 
having the same sign 25 of 29 times (86%). The average of the difference between the 
daily means is 4.0 x 107-°N/m?’, close to the mean values, but well under the standard 
deviation of each. The differences between the Atlas and FNOC curl values are of 
significant size relative to the mean values of each and show considerable scatter, with 
a tendency for the Atlas values to be slightly higher than the FNOC values. However, 
the large variation brings into question the significance of the larger magnitudes of the 
Atlas wind stress curl] values, both for the smaller areas and averaged over one-day 


periods. 
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6. Summary and Conclusions 


There is a demonstrated need for better wind field information over the open 
ocean in order to drive the increasingly complex ocean circulation models. Most 
modelers have relied on fairly simple, artificial wind fields or climatological mean 
fields to provide wind forcing, but the advent of more sophisticated models requires 
realistic wind fields to provide a more highly resolved forcing function. Conventional 
wind data generally lack the combination of resolution and coverage needed for study 
of the temporal and spatial variability of the global ocean. Climatological fields pro- 
vide global coverage, but do not have sufficient spatial or temporal resolution, while 
ship and buoy reports from field investigations provide highly resolved wind data in 
both space and time, but only over limited areas and for short periods of time. The 
conventional fields of most potential use are the synoptic analyses that combine winds 
derived from surface pressure fields and from available ship and buoy reports. Such 
synoptic data is available on a regular, frequent basis over much of the global ocean, 
but is inherently smoothed on the larger, pressure field scales, losing a considerable 
amount of the spatial variability in the wind field. Additionally, in areas such as the 
Southern Ocean, with limited surface pressure and wind reports, even the smoothed 
analyses are questionable. The requirement for a more highly resolved, global spec- 
ification of the surface wind field led to the development of scatterometer theory, as 
a means of inferring surface wind information from a remotely sensed measurement 
of surface roughness from reflected radiation. This initial development culminated in 
the 1978 deployment of an operational scatterometer, SASS, on the SEASAT oceano- 
graphic satellite. Initial evaluations determined that SASS met its stated objectives 
(Jones et al., 1982) and, despite some later qualifications, the three months of data 
returned from this mission provide a unique source of useful scatterometer wind in- 
formation. 

With the deployment of an operational scatterometer, the primary sources of 
more extensive wind information will be the conventional, pressure-derived wind fields 


and scatterometer-derived wind fields. Both wind fields have the potential to provide 
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considerably more information than climatological mean fields, with greatly increased 
temporal and spatial resolution. There are drawbacks and limitations to each, par- 
ticularly in the smoothing of conventional winds on pressure-field scales and in the 
irregular sampling of the scatterometer measurements. However, each is reasonably 
accurate and the question is to select one as forcing for ocean models, or to deter- 
mine the differences, if any, between them. Though the choice of the “best” field in 
an absolute sense is severely hindered by the virtual absence of accurate, widespread 
“surface truth” information, it is as important to evaluate the differences between the 
two fields. An investigation of these differences has been the focus of this study, using 
SASS winds available from R. Atlas at the Goddard Space Flight Center (through 
JPL), and 6-hourly synoptic analysis winds from the Fleet Numerical Oceanographic 
Center (FNOC). 

A one-month subset of the Atlas data, from 12 August to 9 September 1978, over 
the western North Atlantic, from 20° — 50°N and 40° — 80°W, was chosen for this 
study. The comparison required a regridding of the irregularly spaced SASS winds to 
a regular latitude-longitude grid to match the FNOC vector locations. The regrid- 
ding was accomplished by objective mapping, following the procedure of Bretherton, 
Davis, and Fandry (1976). Comparison of the original and mapped Atlas wind fields, 
graphically and statistically, ensured that the mapping accurately reproduced the 
wind field. Comparison of the Atlas and FNOC wind fields to synoptic weather maps 
demonstrated qualitative agreement with surface weather features and general circu- 
lation patterns. Differences between the two fields were examined quantitatively by 
various means, including calculations of field statistics and boxed mean values, scat- 
ter plots of speed, direction, and standard deviation, statistical descriptions of the 
Atlas-FNOC difference field, and calculations of the wind stress curl values for each 
field. 

The Atlas and FNOC fields are similar, in a broad statistical sense, across the 
entire study period, with wide scatter of the individual points about a general trend 
of agreement. The broad correspondence of the two fields is evident in the similar 


distribution of U and V components, for all individual vectors (Figures 11 and 13) 
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and for those averaged over grid areas (Figures 12 and 14). The distribution of 
individual vectors into direction bins also follows a similar pattern for each, though 
more pronounced in the case of the FNOC winds (Figures 11 and 13). The broad 
agreement and large variation of individual points between the two fields is portrayed 
clearly in the Atlas vs. FNOC scatter plots, for both wind speed and direction. 
The relatively low mean values and high variances of the one-day average AU and 
AV components indicate agreement in an average sense, within the limitations of 
individual field variation. The basic similarity of the Atlas and FNOC fields is also 
indicated by comparisons of the wind stress curl fields. The values of the overall 
average curl 7 are very close, the curl 7 values from each field tend to agree in sign 
Over various averaging periods, and the difference between values from each, averaged 
over grid areas as well as one-day periods, is less than the variation of the fields within 
these areas and times. 

Despite this general correspondence, relatively small but consistent differences do 
exist, reflecting the greater smoothing of the FNOC winds on larger scales, and the 
slightly greater Atlas wind magnitudes. The smoothing of the FNOC field is evident 
in the variances of the U and V components, averaged over grid areas, as the the Atlas 
variances are more frequently larger than the FNOC values (Figures 12 and 14). The 
Atlas winds also exhibit a greater range of directions for more of the wind vectors 
(Figures 11 and 13), indicating greater variability in direction. Scatter plots of the 
vector standard deviations (Figure 48) show the greater variance of the Atlas winds 
over one-day and three-day periods, as do the plots of standard deviations of the U 
and V components over the same periods, as well as grouped by latitude (Figures 52 
to 58). The tendency for the Atlas wind magnitudes to be larger shows up in the 
greater numbers of the U and V components (from individual vectors and grid area 
averages) at higher wind speed ranges for the Atlas winds (Figures 11 to 14). Scatter 
plots of all wind speeds from single days, and the scatter plots of all one-day and 
three-day mean wind speeds show that Atlas values are greater than FNOC values 
most of the time. The scatter tilts toward larger Atlas values in 23 of the 29 daily 
plots, 25 of 29 of the mean daily Atlas wind speeds are higher, and all of the three-day 
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mean Atlas wind speeds are higher. The consistency of this trend is striking, but can 
be assigned only limited significance, since in many cases the difference is less than 
the variation of the individual fields over the particular averaging period. Difference 
field statistics show that significant differences exist between individual vectors, but 
that the average differences are small, especially relative to the variation across a 
given averaging period or area. The nearly even division of AU and AV into positive 
and negative values and the relatively even distribution of the difference vector in 
various direction bins (Figure 59) indicates that there is no preferential direction to 
the individual difference vectors. For one-day averages, AU is still split relatively 
evenly into positive and negative values, but AV is negative in most cases, indicating 
that the daily mean Atlas winds are more southerly than the FNOC winds. 

The general similarity and relatively minor differences between the Atlas and 
FNOC fields summarized above characterize most of the data in this study, but some 
extreme differences are present in limited areas. Though significantly different in- 
dividual vectors are scattered throughout the data, regions of large and extensive 
differences are associated in particular with strong circulation features such as hur- 
ricanes. The large differences are due in large part to offsets in the location of the 
low pressure centers, since even relatively small changes can cause large differences to 
appear. The influence of strong circulation patterns is obvious in the difference field 
statistics, where the largest values, for individual points as well as various averaged 
quantities, are found in the vicinity of intense weather, especially Hurricane Ella. The 
mean speed and vector variances of the difference field vectors increase towards the 
end of the study period and also from south to north on average. Both the later 
time (September) and northern regions (40° — 50°N) are associated with relatively 
stronger weather patterns. In addition to these vector differences, the response at 
the highest wind speed range is also significantly different for the two fields, as the 
FNOC values are generally higher. The relatively higher FNOC winds may reflect 
underestimation by the SASS, but the small numbers of vectors involved prevent any 


general conclusion. 
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For this one month period, the Atlas and FNOC winds are similar in their basic 
circulation patterns, in the statistical descriptions of each field, and in the average 
agreement of individual vectors. Most of the differences that do exist lie within the 
considerable scatter due to the variation of each field. However, the extreme differ- 
ences associated with intense circulation features have the potential to greatly reduce 
the overall agreement if their occurrence is frequent enough. The ultimate impact 
on ocean model forcing will depend on the particular area studied and the averaging 
period. In regions of frequent strong weather, the scatterometer and conventional 
winds could yield quite different results, especially at shorter time and space scales. 
Differences at the shorter scales would apply particularly to prediction models for 
forecasting, that depend strongly on the temporal and spatial variations of the in- 
put data. Assuming that the possible underestimation of the highest wind speeds by 
the SASS is resolved in subsequent scatterometers, and that the circulation patterns 
and inferred positioning of weather systems are reproduced more accurately by the 
scatterometer winds (Duffy and Atlas, 1986), the scatterometer information will be 
a significant improvement in these applications. For more quiescent regions and for 
data averaged over longer periods the differences in the two fields should be rela- 
tively minor. The slightly greater spatial resolution of the scatterometer winds will 
be mostly lost in averages over longer periods of time, and should give results similar 


to those of conventional field in those cases. 


64 





Appendix 1: Aspects of atmospheric boundary layer theory 


Motions in the turbulent, atmospheric boundary layer are described by Monin- 
Obukhov similarity theory, which quantifies the balance between turbulence produced 


by wind shear and by buoyancy. A different velocity scale is important in this theory: 
: reer 
us = (7/9)? =| uw |? 


OT, 


T = p(us)* 


where: 


ux = friction velocity (m/s) 


The friction velocity, u,, is a surface flux parameter used as another means of repre- 
senting wind stress, in this case the kinematic stress. Also introduced is the Monin- 
Obukhov scaling length, L, which gives the height at which the turbulent energy 
production from non-neutral buoyancy and from mechanical shear in the neutral case 
are equivalent (Large and Pond, 1981). Essentially, it is the height at which shear 
and buoyancy effects are approximately the same. L can be expressed as: 


Ted epen mn (Cri gla lic 
egH ~~ «g(Ts — Ta) 


where: 
T, = air temperature at surface (K) 
Cp = Specific heat capacity of air (1030 J/kgK) 
« = von Karman’s constant (~ .35) 
g = acceleration of gravity (9.8m/s”) 
H = flux of sensible heat (W/m?) 
Uio = mean wind speed at 10m above mean sea level (m/s) 
; = mean sea surface temperature (K) 


(Stewart, 1985) 
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In both of these forms it is assumed that temperature effects are much greater than 
those of humidity, which is generally true for mid-latitudes, but does not generally 
hold in the tropics or under conditions of high evaporation, requiring adjustment to 
the equations in those cases (Stewart, 1985). A more general form can be found in 
Pond et al. (1971). For typical values of air-sea temperature difference, L is usually 
greater than 30 meters, with the possible exception of very light wind conditions. 
Since this height is above the area most crucial to scatterometry, stability plays a 
secondary role and enters as adjustments to profiles derived by assuming a neutrally 
stable atmosphere. 

Neglecting stability, the velocity profile is logarithmic with height. The stability 
corrections are functions of the non-dimensional stability parameter, 2 , where z is 


L 
the height above the surface and: 


2 — 0 neutral stability 
é >0 stable conditions 
2 <0 unstable conditions 


This yields the mean wind profile, in gradient form: 


Kz du 


on = > 


uy dz 


where: 
z = height above mean sea level (m) 
Om = om(#) = dimensionless wind shear 


(Businger et al., 1971) 


Integrating this yields the actual profile: 


itl Z Z 
Un = (n(2) + bm(2)) 
where: 
U, = wind speed at some height z (m/s) 


Z. = roughness parameter (m) 
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Since ¢,,(0) = 1, the neutral wind profile is logarithmic with height, as required. The 


form of ¢,, is known only empirically, with one formulation as: 


bm = (1— 152)-* (¢ 
om = (14+ 4.74) (2 > 0) 
Pm = 1 (= 


(Businger et al., 1971) 


Stewart (1985) shows that if the influence of stability is relatively weak, ¢,, can be 


expanded in a power series, yielding the profile solution: 


ae Z, Z— Zo 
U, = —(In(— sone 
»= San(2) + (45%) 
where: 
B=4.T 2>0 
B= 3.8 2 <0 


The major effect of stability is simply to shift the velocity profile away from the 
logarithmic form at non-neutral conditions. To estimate stress with a drag law from 
this profile, either the velocity or drag coefficient must be corrected for stability. 
Given a similar equation for the virtual temperature profile and the expression 
for L, we have a closed set of equations that can be solved for the fluxes given the 
profiles, or vice versa. However, the solution depends on finding a suitable expression 
for the roughness parameter, z,. (Halberstam, 1980). This parameter arises from the 
mathematical process and can be considered a statistical description of the density 
and height of the roughness elements (Wu, 1980). In the equations, z,. represents 
the height, close to the surface, at which U vanishes. In reality the velocity does not 
vanish, because, in the few centimeters immediately above the surface, viscosity effects 
become important, and the above relations are no longer valid. For solid surfaces z, 
is closely related to surface roughness and usually taken as constant, but a constant 


Z. is not applicable to the moving ocean surface. Since wave heights (which form 
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the roughness elements) are changed by the wind, z. has a pronounced wind speed 
dependence. The exact form of this dependence is not known, and many studies 
have been done to specify an empirical relation between z, and u,. Charnock (1955) 
assumes that the shorter, steeper waves are primarily responsible for variations in the 
relevant roughness parameter and that these are usually in equilibrium with the local 


wind. His predicted form for z, was: 


mi 





Di 


where: 


a = a proportionality constant (0.0156) 


This information fits the observations well and also fixes the form of the drag coeffi- 
cient, since Cp is related to u,. The relation is rationalized physically by Wu (1968), 
as an equation of state that characterizes the equilibrium between wind and waves 
with gravity waves acting as roughness elements (Wu, 1980). Further refinements to 
this form have been made, both those that assume a unique relation (Cardone, 1969; 
Garratt, 1977) and those that also incorporate further adjustments for swell and wind 
fetch (Kitaigorodskii and Zaslavskii, 1974). Garratt (1977) reviews several of these 


and concludes that a simple z, -to-u, relation is adequate, of the form: 
Zee calm ie (CGS units) 


Halberstam (1980) showed that the u* derived from different forms of z, did not 
correlate significantly differently with backscatter measurements. The absence of 
significant differences is mostly a function of the considerable scatter in the SASS 
backscatter data, since it tends to swamp the differences due to different z, formula- 
tions. 

As mentioned above, this parameterization of z, also fixes the form of the drag 


coefficient, since Cp = cael then using the derived form of z, in the log profile: 


K 
Cp — ines)! (Stewart, 1985) 


3 
aus 
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With known values of kK, g, z, and a, these equations can be solved numerically to 
yield Cp(Uio) (Stewart, 1985). Many measurements have been made to determine 
an empirical Cp-to-U,.; relation, both through direct and indirect techniques (ie: 
Smith, 1980; Large and Pond, 1981). The specification of wind speed dependence has 
been a major issue. Measurements made by Smith (1980), Large and Pond (1981), 
and others show a definite, though widely scattered, increase of the drag coefficient 
with wind speed. Wu (1980) uses scaling laws for Cp to explicitly demonstrate that 
Cp must increase with wind speed if both the wind speed follows a logarithmic profile 
and the Charnock relation holds. Physically, this increase, especially at higher wind 
speeds, can be tied to the different character of the ocean surface, as compared to 
a solid surface. Over a solid surface, as wind speed increases the viscous sublayer 
becomes thinner and roughness elements begin to protrude through it, increasing the 
roughness length (and drag coefficient). In the fully rough condition, the elements 
are completely exposed and since the roughness length will not increase with further 
increase of velocity, the drag coefficient is constant. On the other hand, on the 
ocean surface the short waves are the roughness elements and these continue to grow 
with increasing wind speed, even though a fully rough condition exists. So the drag 
coefficient continues to increase with increasing wind velocity. One notable aspect of 
ocean surface roughness is the presence of two markedly different regimes for lower 
and higher wind speeds. Within the atmospheric boundary layer, under neutral or 
unstable stability conditions, instabilities and large-scale coherent structures develop 
for winds higher than about 7 m/s (Brown, 1986). This change of behavior in the 
boundary layer flow is reflected by a discontinuity in plots of surface roughness versus 
backscatter. This pattern has been recognized for quite some time - see Woiceshyn et 
al. (1986) for a more complete discussion. 

The discontinuous behavior of Cp cannot be fit with a single linear relation, in 
other words, by choosing a constant drag coefficient over the entire range of wind 
speeds. The alternative generally chosen is to fit two linear relations, across lower 
and higher wind speed regimes. Large and Pond (1981) have shown that a form of 


the drag coefficient that is constant at low wind speeds and varies linearly with U at 
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higher wind speeds fits the observations best. They relate the neutral drag coefficient 
referenced to 10 m (Cpnio) and the wind speed at 10 m (Uj): 


lee 10> Uio < 11m/s 


Cpn10 = (G40) 065Uiq) < 10m: Wig > iim/s 


They also provide a method for calculating stress from wind speed at other heights, 
using this formula and the air-surface temperature difference. The variability of the 
drag coefficient with wind speed affects transient forcing as well as the mean wind 
stress. Variable drag coefficients have been shown to yield transient forcing that is 
approximately 30% higher than that for constant coefficients (Smith, 1980). The 
increased contribution of variable drag coefficients is an important consideration in 
selecting the forcing parameterizations and U-to-7 conversion, if necessary, for ocean 
modeling. With these relations, wind speed at various levels in the surface layer can 
be related to wind speed near the surface, and to stress on the surface through a drag 
coefficient. Stewart (1985) outlines the procedures: 1) Given the wind velocity at 
some height, U,, and the air-sea temperature difference, find L and estimate é Zee 
> is negligible and stability effects can therefore be ignored, use the logarithmic wind 
profile and Charnock’s value of z, to find U;¢¢ (usually at 10 m); 3) If Z is small, but 
not negligible, correct U,, to the value it would have if = = 0 (at the same height), 
then find U;e¢ with the log profile; 4) If é is not small, then U;er can be estimated 
by numerical integration of the velocity profile, but it is generally better to get a new 
measurement of U, at a lower height. The extrapolations made for large é values are 


not very accurate. 
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Appendix 2: Ocean waves and energy transfer 


Capillary waves have wavelengths of about two centimeters, at which scale sur- 
face tension is a dominant influence. The sea surface can be described with spectra 
calculated from a three-dimensional Fourier transform, and the spectrum of wind- 
generated gravity waves is well known (Stewart, 1985). Spectra for the higher fre- 
quency capillary waves are not as well understood, nor is the actual spectra of the 
entire sea surface. Lleonart and Blackman (1980) discuss spectral characteristics of 
capillary waves and summarize some high-frequency wave measurements. Their stud- 
ies indicate that the shape of the frequency spectrum, over the capillary wave range, 
depends on friction velocity, viscosity, and surface tension, and therefore on the local 
wind field. The height of small wavelength waves depends on balancing the input 
of energy from the wind and from other waves through wave-wave interactions with 
the loss of energy through viscosity, surface films, and other wave-wave interactions 
(Stewart, 1985). The method by which wind energy is input to capillary waves is still 
not completely understood, as discussed by Phillips (1977). Variations in both sur- 
face pressure, a normal force, and shear stress, a tangential force, contribute to wave 
growth. Phillips (1977) makes a distinction between two types of stress fluctuations, 
those generated by atmospheric turbulence and those induced by flow over the wavy 
surface. Turbulent stresses, due to their random nature, tend to contribute energy 
across a large frequency range, while the induced stresses contribute energy more se- 
lectively, enhancing the growth of only particular wave components (Phillips, 1977). 
An important issue is the allocation of the total momentum transfer, between the 
fraction associated with tangential stress, which is assumed to correlate with capil- 
lary waves, and the portion associated with the normal force of form drag, acting on 
the longer wave components (Stewart, 1985). Kitaigorodskii (1970) summarizes some 
of this work and finds that 70-85% of the momentum transfer to waves is linked to 
capillary waves, and the rest to longer waves. Theoretical calculations (Brooke Ben- 


jamin, 1959; Miles, 1962) and experimental observations (Kendall, 1970) indicate that 
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the normal stresses are more important than the tangential stresses in energy trans- 
fer. Pierson et al. (1986) cite this work in their choice of a new physical parameter to 
relate to microwave backscatter, discussed in Section 2. The presence of significant 
energy transfer to the longer waves means that this long wave energy input must re- 
late nearly constantly to the short wave energy input in order for the assumed linear 
wind-to-short wave spectra relations to hold (Brown, 1986). While an approximately 
constant relation holds for the near steady state conditions commonly found over the 
ocean, it does not hold under non-steady state conditions, such as those found in the 
vicinity of fronts and other weather systems. Phillips (1977) formulates an expression 
for the air-ocean energy flux that portion energy into wave motion and currents. It 
includes terms for the flux of energy into waves alone, for energy flux from normal and 
tangential stress variations, for the loss of mean flow energy from Reynold’s stresses 
of the induced motion working against the mean velocity gradient, and for energy 
loss from wave-induced motion from molecular viscosity or variations in the turbulent 
Reynold’s stress (Phillips, 1977). 

In addition to the role of long waves in the partition of transferred energy, the 
presence of these waves may attenuate or enhance short wave growth, depending 
on the relative magnitude and orientation of the two (Phillips, 1978). The short 
wave shape, and therefore the associated surface roughness, can also change due to 
wave-current interactions, under the same wind conditions (Longuet-Higgins, 1978). 
Rain can significantly alter the short wave spectra, primarily by causing damping at 
higher rates of fall. Sea surface temperature variations affect short waves through 
the associated changes in viscosity, up to a factor of two across a 0-30°C temperature 
range (Pierson et al., 1986). Mixed layer dynamics, such as those responsible for 
Langmuir circulations can also significantly affect short wave spectra (Brown, 1986). 
These effects illustrate the complex system of multiple, interacting processes that 
complicate the basic picture of energy flux to the waves as a tangential stress correlated 
with short wave elevation and a normal stress in phase with the slopes of the longer 


Waves. 
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Appendix 3: Aspects of radar scatter theory 


Microwave scatterometry is based to a large extent on known principles of radar 
scatter from rough surfaces. Calculating the backscatter from incident radiation on 
rough surfaces requires matching the electric and magnetic fields of the radiation 
across a boundary of known shape. Specifically, the scattering cross-section, o, can 
be written as 4m times the ratio of the scattered radiant intensity at the receiver to 
the power density incident on the area: 

o = 4nAcos 9(25) 

0; 
where: 
o = backscatter cross-section (m7) 
A = surface area (m7) 
@ = incidence angle measured from vertical 
I, = radiant intensity (W/sr) 
OQ; = radiant flux (W) 
(Stewart, 1985) 


The incidence angle, @ , is the angle between the incident radar beam and the local 
surface vertical. Power density is the radiant flux, O;, divided by the projected area, 
Acos(@). Using the radar equation, an expression for the received power can be found 
in terms of the transmitted power and the scattering cross section. Stewart (1985) 
develops an expression of this type for a point target. Since the sea is a distributed, 
rather than a point target, the chosen radiation measure is the normalized radar 
backscatter cross section (NRCS), o° . The dimensionless 0° is a ratio of the reflected 
to incident energy across a unit surface area, rather than the projected area. From 
these relations we can, in theory, find the power scattered toward the radar receiver 
if the form of o° is known. So the problem is to find o° for a given sea surface. 
Current approaches to wave scattering from statistically rough surfaces fall into 


two categories: 1) specular or tangent-plane theory for slightly wavy surfaces, and 
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2) Bragg-resonance theory for slightly rough surfaces. These two physical mecha- 
nisms take place over different ranges of incidence angles, with specular scatter at 
near vertical 8, and Bragg scatter at larger angles. Specular scatter is a mirror-like 
reflection from wave facets that are oriented to reflect energy directly back toward the 
radar receiver. Because the direct reflection comes from wave slopes perpendicular 
to the incidence angles, and the fact that ocean wave slopes are rarely more than 
20°, this mechanism is restricted to a small range of incidence angles. To apply this 
mechanism, wave facet shape must be specified, an estimate made of the probability 
of proper facet orientation made, and then the energy reflected by the facet calculated 
(Stewart, 1985). The estimate is generally made by use of a tangent plane approxima- 
tion, through either a physical optics or geometrical optics approach. Stewart (1985) 
discusses the process in more detail, with expressions that sum the scatter over the 
incidence area, weighted by a probability distribution of surface slopes. The sea sur- 
face slopes are assumed to have a Gaussian distribution that is anisotropic about the 
wind direction. Those ocean wave lengths that are shorter than the radar wavelength 
do not contribute to this type of scatter. The return from specular scatter is isotropic 
in azimuth, which means that only magnitude and not direction information can be 
determined. 

At larger incidence angles, the number of facets available for direct scatter de- 
creases rapidly and another mechanism dominates the scatter. This is Bragg scatter, 
a resonant reflection from those wave components that match the projected radar 
wavelength on the ocean surface. It can be described as the combination of scattering 
from different surface elements, which, through constructive interference enhances the 
scattering from regular surface structures with favorable wavelengths, and through de- 
structive interference diminishes all other reflections. The physical principle behind 
this type of scatter is that of a diffraction grating, and the radiation scatters in di- 
rections set by spacing and orientation of the regular surface elements relative to the 
radar wavelength (Stewart, 1985). The ocean surface can be viewed as a superposition 
of plane waves, and considering scattering as a linear process, then each individual 


component acts as a diffraction grating for a particular radiation wavelength. Since 
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this resonant scatter adds coherently in the far field, with power given by Bragg’s 
equation, it is known as Bragg scatter. Resonant scatter theory can be used with a 
physical optics approach to relate the backscatter measurement to the ocean surface. 
In a formulation developed by Wright (1968), using linearized boundary conditions, 
the backscatter cross section is directly proportional to the mean spectral density of 


the short-wavelength ocean waves: 
o°(9,k)3; = 167k? cos* 8 | g;;(0,€,) |? U(kx, ky) 


where: 
W(kx, ky) = 2-D Fourier spectrum of surface height fluctua- 
tions 
gij(9,€,) = first-order scattering coefficients 
€, = complex dielectric constant of seawater 
k, = radar wavenumber 


6 = incidence angle 


The actual dependence is on the spectrum of ocean wave slopes, but since that is 
linearly related to the surface height spectrum (ie: to the wave height spectrum), one 
can be used in place of the other (Stewart, 1985). This form combines the spectral 
densities of waves traveling toward and away from the antenna, since their respective 
doppler shifts are indistinguishable from each other. 

For resonant scatter theory to hold, one condition is that k,h < 1, where h is 
the wave height. This condition is not satisfied for the short radar wavelengths used 
with capillary waves, but the scatter from these waves is handled by considering the 
ocean as a composite surface (Stewart, 1985). Composite surface theory assumes that 
the ocean surface has two separable scales of motion, with short-wavelength ripples 
superimposed on the longer gravity waves (Figure 66). Wright (1968) and Bass et 
al. (1968) first developed the basic principles of composite surface scattering theory, 
characterizing backscatter from the actual sea surface over the Bragg scatter range of 


incidence angles. Over small areas the Bragg scatter is calculated from the local wave 
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field and local surface orientation. A probability distribution based on the slopes of the 
longer waves can then be used to integrate the local scatter values over a larger area. 
These mechanisms valid at intermediate incidence angles are then combined with the 
specular scatter mechanism at small incidence angles to produce a backscatter model 
for the ocean over a large range of incidence angles (Valenzuela, 1978). Laboratory 
and field experiments (see Stewart, 1985 for partial summary) support the two-scale 
Bragg scattering model over incidence angle ranges of approximately 20° < @ < 70° 
for vertical polarization and 20° < @ < 60° for horizontal polarization. For 6 < 20°, 
specular scatter starts to interfere, and it becomes the dominant scattering process for 
incidence angles on the order of the rms wave slope, roughly 5°-8° for wind generated 
waves. In the 10°-20° incidence angle range, the backscatter is relatively insensitive 
to wind speed. At the largest incidence angles (9 > 60°), wedge diffraction may 
become an important mechanism (Lyzenga et al., 1983), and shadowing by larger 
waves may prevent signal reception. Also, at these large angles, the intermittent signal 
scattered from white caps may be as strong as that from Bragg scatter, significantly 


contaminating the return. 
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Table 1: Ungridded Atlas SASS component wind statistics 
Individual area averages 
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DAY 


224 U 


225 U 


226 U 


Ze 


228 U 


Ot 


230 U 


231 U 


232 U 


233 U 


234 U 


235 U 


236 U 


237 U 


238 U 


239 U 


240 U 


241 U 


242 U 





-0.720 
2.236 
-1.897 
1.049 
-2.776 
-0.204 
0.611 
-0.789 
Zee? 
-1.102 
-0.322 
0.809 
-0.375 
1.828 
3.052 
2.170 
0.113 
1.037 
0.012 
-1.101 
-1.103 
-1.076 
0.316 
0.925 
1.680 
0.936 
0.970 
-1.369 
-1.435 
2.623 
-1.515 
2.966 
-2.549 
1.997 
-3.267 
0.178 
-2.858 
-2.144 


VAR 


(m?/s*) 


26.060 
10.957 
32.461 
10.627 
32.242 
12.890 
45.880 
16.347 
45.433 
16.523 
29.390 
10.624 
23.460 
10.026 
6.765 

told 

32.800 
19.621 
30.609 
25.519 
25.662 
23.019 
29.565 
16.703 
27.300 
15.692 
26.786 
7.525 

33.986 
Ze ES 
21.785 
14.645 
14.509 
10.473 
15.535 
17.292 
17.301 
17.253 
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WVAR 





0.982 
0.688 
0.902 
0.908 
0.808 
0.998 
0.994 
0.965 
0.912 
0.933 
0.998 
0.943 
0.995 
0.751 
0.422 
0.617 
1.001 
0.949 
1.001 
0.956 
0.956 
0.953 
0.998 
0.952 
0.907 
0.948 
0.969 
0.803 
0.944 
0.755 
0.907 
0.626 
0.692 
0.725 
0.593 
1.000 
0.681 
0.791 





5.105 
3.310 
5.697 
3.260 
5.678 
3.590 
6.773 
4.043 
6.740 
4.065 
5.421 
3.259 
4.844 
3.166 
2.601 
2.742 
Set 26 
4.430 
5.533 
5.052 
5.066 
4.798 
5.437 
4.087 
9.225 
3.961 
5.176 
2.743 
5.830 
4.595 
4.667 
3.827 
3.809 
3.236 
3.941 
4.158 
4.159 
4.154 


Table 2: Ungridded Atlas SASS component wind statistics - 1-day averages 
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250 U 


251 U 


252 U 


-2.497 
-2.421 
-2.476 
-0.019 
-2.399 
-1.218 
-1.492 
0.741 
-0.710 
2.888 
0.805 
4.175 
-1.746 
1.782 
-1.885 
-0.397 
1.109 
3.620 
0.184 
2.454 


15.265 
12.829 
12.786 
24.092 
29.893 
11.975 
39.713 
29.367 
31.102 
34.086 
41.735 
28.833 
40.097 
26.404 
34.931 
18.012 
29.008 
53.994 
51.466 
62.263 
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0.713 
0.689 
0.677 
1.003 
0.820 
0.892 
0.949 
0.984 
0.986 
0.805 
0.986 
0.624 
0.931 
0.894 
0.911 
0.995 
0.961 
0.806 
1.001 
0.914 


3.907 
3.582 
3.576 
4.908 
5.089 
3.461 
6.302 
9.419 
5.577 
5.838 
6.460 
5.370 
6.332 
9.138 
9.910 
4.244 
5.386 
7.348 
7.174 
7.891 


169 
169 
396 
396 
371 
371 
456 
456 
523 
523 
765 
765 
905 
905 
266 
266 
510 
910 
474 
474 





U component 


REV. 





658 Ul 
U2 
664 U1 
U2 
Us 
665 U1 
672 Ul 
U2 
U3 
678 U1 
679 Ul 
U2 
686 U1 
WZ 
692 U1 
693 U1 
wi 
Ls 
700 U1 
701 U1 
U2 
U3 
708 U1 
U2 
715 U1 
WZ 
U3 
U4 
716 Ul 
721 Ul 
2 
2 
729 U1 
U2 
735 Ul 





4.098 
-4.553 
-7.179 
-4.128 
5.628 
-0.152 
6.097 
4.340 
-5.128 
-1.757 
-7.142 
-1.396 
-3.526 
-8.387 
-2.764 
-7.849 
-0.125 
7.091 
-4.162 
9.221 
2.582 
-6.729 
-4.737 
4.551 
12.911 
8.060 
2.980 
-4.430 
3.059 
6.785 
-5.810 
4.977 
2.221 
-6.259 
10.185 


VAR 


(m*/s*) 


3.378 
2.720 
3.022 
0.864 
14.099 
3.945 
4.069 
2.850 
0.858 
2.174 
1.662 
Zl 
1.159 
1.265 
6.509 
4.805 
4.279 
2.947 
2.910 
1.918 
3.948 
9.848 
7.618 
2.404 
6.562 
i250 
9.355 
19.636 
3.924 
8.259 
11.100 
8.061 
7.029 
S.4 2 
10.020 
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0.168 
0.116 
0.055 
0.048 
0.311 
1.022 
0.099 
0.132 
0.032 
0.417 
0.032 
0.942 
0.085 
0.018 
0.466 
0.072 
1.025 
0.055 
0.144 
0.066 
0.376 
0.115 
0.255 
0.104 
0.040 
0.019 
0.594 
0.505 
0.276 
0.152 
0.249 
0.279 
0.595 
0.088 
0.088 





1.838 
1.649 
1.738 
0.930 
3.755 
1.986 
2 Oey, 
1.688 
0.926 
1.475 
1.289 
4.666 
OTs 
1.125 
Zool 
2.192 
2.069 
Lilt 
1.706 
1.385 
1.987 
2.418 
2.760 
1.551 
2.962 
i2t 
3.059 
4.431 
Sei 
2.874 
3.332 
2.839 
2.651 
1.942 
3.165 


Table 3: Mapped Atlas SASS component wind statistics - area averages 
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U3 
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U3 
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O71 22 
1.196 
2.048 
-4.737 
-3.059 
-1.584 
2.441 
2.106 
3.608 
4.288 
2.241 
-5.407 
1.473 
7.407 
-5.084 
3.120 
-1.961 
5.276 
-4.670 
5.699 
-1.857 
-5.387 
6.374 
3.760 
-2.134 
3.461 
-2.618 
3.985 
0.283 
-0.575 
-3.318 
2.142 
-6.572 
-5.637 
-2.672 
1.621 
eeiah 
-4.865 
267 
6.730 
-0.472 
1.281 


6.327 
0.995 
9.151 
17.283 
1.235 
10.466 
18.362 
6.991 
fe 
18.090 
1.346 
1.577 
0.579 
PO 
16.807 
2.350 
5.975 
3.874 
45.232 
5.8953 
8.206 
11.789 
4.105 
1.843 
7.134 
5.845 
40.876 
1.860 
2.658 
43.080 
17.538 
eet 
5.832 
7.018 
5.654 
1.204 
30.229 
16.976 
7.704 
2.041 
6.441 
eo ly, 
14.181 
33.791 
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0.308 
0.200 
1.018 
0.931 
0.160 
0.319 
0.666 
0.747 
0.552 
0.817 
0.094 
0.079 
0.104 
0.043 
0.902 
0.041 
0.189 
0.286 
0.946 
0.174 
0.275 
0.268 
0.552 
0.060 
0.150 
0.295 
0.923 
0.135 
0.281 
0.742 
1.016 
1.012 
0.349 
0.615 
0.116 
0.037 
0.822 
0.878 
0.945 
0.080 
0.814 
G22) 
1.012 
0.961 


2.515 
0.997 
2.270 
4.157 
Pit 
3.235 
4.285 
2.644 
2.685 
4.253 
1.160 
1.256 
0.761 
1.141 
4.100 
1.533 
2.444 
1.968 
6.725 
2.419 
2.865 
3.434 
2.026 
1.357 
2: 6ucl 
2.418 
6.393 
1.364 
1.630 
6.564 
4.188 
4.614 
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2.649 
2.378 
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2eishO 
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9.813 


49 
49 
49 
121 
121 
121 
aa 
49 
64 
49 
49 
49 
36 
49 
49 
49 
36 
64 
36 
49 
49 
49 
36 
36 
49 
36 
36 
49 
49 
49 
49 
36 
49 
36 
36 
36 
49 
64 
49 
49 
49 
64 
36 
PAI 





'? 


U2 
U3 
836 U1 
U2 
837 Ul 
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U2 
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U2 
U3 
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866 U1 
873 U1 
U2 
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U2 
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901 U1 
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U2 
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U2 
U3 
930 U1 
U2 
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WZ 
U3 
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U2 
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U2 
U3 
973 U1 
U2 
U3 
987 Ul 
988 Ul 
we2 
1002 U1 
U2 
U3 


4.800 
1.204 
5.337 
0.872 
-1.211 
4.956 
-2.003 
2.251 
-5.150 
4.101 
-4.587 
0.019 
3.035 
-4.072 
6.919 
-0.940 
2.064 
-2.878 
-3.852 
-3.899 
-1.488 
-3.649 
2.3907 
-6.806 
-3.508 
0.467 
-5.344 
-3.284 
1.080 
2.141 
-2.710 
-2.246 
-2.170 
-2.754 
-5.150 
3.063 
-1.423 
-6.993 
-5.437 
6.973 
-1.070 
4.029 
5.439 
-3.103 


7.074 
2.843 
9.696 
45.911 
2.185 
1.180 
17.167 
5.439 
6.097 
4.562 
29.789 
36.865 
9.204 
19.614 
5.432 
14.347 
7.602 
8.729 
292 
5.717 
2.655 
11.407 
11.493 
2.676 
6.299 
6.682 
2.066 
20.237 
1.115 
0.602 
1.366 
10.658 
0.891 
0.297 
4.594 
4.395 
0.671 
4.408 
6.194 
17.335 
L6ers 1 
2.348 
6.236 
7.825 
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0.236 
0.669 
0.256 
0.992 
0.604 
0.046 
0.824 
0.513 
0.187 
0.214 
0.590 
1.029 
0.366 
0.546 
0.102 
0.949 
0.544 
0.518 
0.076 
0.274 
0.551 
0.463 
0.681 
0.055 
0.341 
0.988 
0.068 
0.664 
0.495 
0.116 
0.157 
0.688 
0.160 
0.038 
0.148 
0.322 
0.251 
0.083 
0.173 
0.265 
0.945 
Ome i 
0.175 
0.450 


2.660 
1.686 
3.114 
6.776 
1.478 
1.086 
4.143 
2.332 
2.469 
2.136 
9.458 
6.072 
Peco 
4.429 
2.331 
3.788 
Ziebart 
2.955 
1.101 
2.391 
1.630 
S.017 
3.390 
1.636 
2.910 
2.085 
1.437 
4.499 
1.056 
0.776 
1.169 
3.265 
0.944 
0.545 
2.143 
2.096 
0.819 
2.100 
2.489 
4.164 
4.097 
1.532 
2.497 
De OW 


Sil 
64 
36 
121 
64 
36 
49 
121 
64 
64 
81 
36 
49 
81 
49 
121 
36 
49 
36 
121 
49 
121 
64 
49 
49 
49 
49 
36 
36 
36 
49 
49 
49 
36 
49 
36 
36 
49 
121 
36 
100 
121 
49 
121 





1003 U1 
1008 U1 
U2 
1009 U1 
U2 
1016 U1 
U2 
1017 U1 
1030 U1 
1031 U1 
WZ 
U3 
1051 U1 
1052 U1 
U2 
1059 U1 
U2 
1060 U1 
1074 U1 
U2 
U3 
U4 


0.150 
3.190 
6.538 
-1.800 
4.763 
3.393 
-7.608 
2.008 
-4.464 
5.608 
1.669 
-5.088 
3.665 
1.631 
5.786 
0.816 
-2.113 
-1.723 
-12.984 
6.019 
5.368 
-1.800 


V component 


REV 


658 V1 
V2 
664 V1 
V2 
3 
665 V1 
672 V1 
V2 
V3 
678 V1 
679 V1 
V2 
686 V1 
V2 
692 V1 
693 V1 


MEAN 
(m/s) 


2.348 
0.783 
0.697 
0.504 
6.018 
6.552 
1.903 
5.178 
-0.345 
2.290 
0.300 
4.115 
-1.647 
-0.190 
4.238 
-0.022 


3.905 
7.833 
68.106 
10.402 
3.753 
13.884 
9.405 
2.917 
12.801 
4.806 
21.126 
5.881 
46.865 
16.638 
21.046 
16.569 
5.648 
0.750 
16.989 
86.731 
25.262 
6.988 


VAR 
(m*/s*) 


5.671 
0.763 
7.154 
0.670 
9.940 
8.178 
4.385 
1.530 
2.044 
2.599 
1.428 
5.656 
6.412 
4.279 
2.136 
8.535 
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1.010 
0.440 
0.618 
0.767 
0.142 
0.553 
0.140 
0.425 
0.392 
0.133 
0.906 
0.186 
0.782 
0.870 
0.389 
0.981 
0.565 
0.203 
0.092 
0.716 
0.472 
0.697 





0.512 
0.563 
0.961 
0.733 
Os ZK 
0.161 
0.556 
0.054 
0.971 
0.330 
0.966 
0.252 
0.713 
1.012 
0.107 
1.021 


1.976 
2.199 
8.253 
3.225 
1.937 
3.726 
3.067 
1.708 
3.578 
2.192 
4.596 
2.425 
6.846 
4.079 
4.588 
4.070 
2.307 
0.866 
AZZ 
9.313 
5.026 
2.643 


SDEV 


(m/s) 


2.981 
0.873 
2.675 
0.818 
3.153 
2.860 
2.094 
eo) 
1.430 
1.600 
1.195 
2.378 
2.532 
2.069 
1.462 
2.921 


64 
36 
121 
121 
121 
49 
81 
36 
rn 
64 
36 
36 
121 
100 
49 
49 
49 
36 
49 
49 
49 
36 


49 
36 
36 
64 
36 
36 
36 
49 
36 
49 
36 
36 
49 
49 
36 
49 


Oe penanD 





V2 
V3 
700 V1 
701 V1 
2 
V3 
708 V1 
V2 
715 V1 
V2 
V3 
V4 
716 V1 
721 V1 
G2 N 1 
We 
729 V1 
V2 
735 V1 
736 V1 
V2 
V3 
744 V1 
V2 
V3 
750 V1 
V2 
V3 
751 V1 
758 V1 
772 V1 
773 V1 
779 V1 
4 
V3 
780 V1 
787 V1 
V2 
V3 
V4 
793 V1 
V2 
794 V1 
V2 


2.856 
5.321 
-9.105 
-2.044 
3.107 
0.644 
1.106 
0.611 
-2.938 
-0.544 
-0.121 
-1.537 
2.665 
-1.090 
0.365 
-2.375 
-1.376 
-0.686 
-0.944 
-0.725 
-3.226 
4.880 
0.242 
3.527 
1.789 
-0.649 
-0.735 
1.711 
3.950 
3.703 
3.010 
-0.524 
-0.439 
3.342 
0.307 
2.052 
-4.738 
4.819 
2.837 
1.296 
-5.169 
-6.090 
2.040 
0.223 


0.576 
1.389 
1.591 
1.692 
1.875 
1.992 
3.433 
1.833 
33.891 
4.302 
0.544 
10.174 
1.623 
6.569 
4.426 
5.131 
1.944 
3.524 
3.958 
2.996 
1.533 
0.395 
7.246 
4.739 
Zeon 
2.233 
3.176 
9.999 
5.721 
3.839 
5.429 
2.342 
1.081 
2.003 
35.855 
6.889 
28.877 
7.055 
8.013 
5.452 
1.560 
11.234 
6.962 
2.434 
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0.066 
0.047 
0.058 
0.291 
0.163 
0.847 
0.753 
0.850 
0.810 
0.954 
1.001 
0.825 
0.187 
0.858 
0.998 
0.481 
0.512 
0.904 
0.835 
0.846 
0.129 
0.016 
1.000 
0.276 
0.412 
0.847 
0.870 
0.783 
0.270 
0.220 
0.378 
0.918 
0.864 
0.184 
1.018 
0.632 
0.568 
0.235 
0.504 
0.777 
0.055 
0.234 
0.527 
1.000 


0.759 
1.179 
1.261 
1.301 
1.369 
i 
1.853 
1.354 
5.822 
2.074 
0.737 
3.190 
1.274 
2.963 
2.104 
2.265 
1.394 
1.877 
1.990 
1.611 
1.238 
0.628 
2.692 
Zale 
1.492 
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Table 4: Mapped Atlas SASS component wind statistics - 1-day averages 
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DAY MEAN VAR WVAR  SDEV N 
(m/s) (m?/s?) (m/s) 
224 U -0.109 27.102 1.002 5.206 378 
V 2.548 9.800 0.603 3.130 378 
225 U -2.146 26.860 0.856 5.183 340 
V 1.444 8.633 0.807 2.938 340 
226 U -2.580 26.378 0.800 5.136 463 
V -0.201 14.313 0.999 3.783 463 
227 U 1.222 41.364 0.967 6.432 499 
V -1.085 12.714 0.917 3.566 499 
228 U 2.929 41.766 0.831 6.463 453 
V -0.926 9.551 0.919 3.090 453 
229 U 0.077 23.575 1.001 4.855 631 
V 0.938 8.676 0.909 2.946 631 
230 U -0.189 19.515 1.000 4.418 695 
V 1.500 7.567 0.772 2.751 695 
231 U 3.490 1.869 0.133 1.367 134 
V 2.314 7.031 0.570 2.652 134 
232 U 0.937 Psos 0.971 5.249 466 
V 0.928 18.717 0.958 4.326 466 
233 U 1.256 27.875 0.948 5.280 636 
V -0.878 20.890 0.966 4.571 636 
234 U -0.069 22.367 1.001 4.729 623 
V -0.661 16.108 0.975 4.013 623 
235 U 0.930 25.584 0.969 5.058 662 
V 1.383 14.450 0.884 3.801 662 
236 U 1.728 22.131 0.882 4.704 693 
V 1.755 15.901 0.839 3.988 693 
237 U 1.730 12.742 0.813 3.570 206 
V -2.126 4.021 0.472 2.005 206 
238 U -0.722 30.494 0.985 5.522 545 
V 3.213 15.631 0.603 3.954 545 
239 U -0.785 24.319 0.978 4.931 372 
V 3.671 8.444 0.386 2.906 372 
240 U -2.849 10.468 0.564 3.235 412 
V 1.587 8.519 0.773 2.919 412 
241 U -2.960 13.620 0.609 3.690 453 
V -0.269 20.098 0.999 4.483 453 
242 U -2.283 18.757 0.785 4.331 260 
V -2.884 17.582 0.681 4.193 260 
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Table 8: FNOC component wind statistics - area averages 


U Component 


DAY/TIME MEAN VAR WGTVAR SDEV 
(m/s) (m%*2/s%2) (m/s) 

22e°00 Ul Zao 3.364 0.344 1.834 
Wiz -3.856 023.56 O20Z5 02621 

yza 12 Ul -/.366 02} Ors Le O10 
WZ -2.861 Zeooe 0.260 1 62 

U3 2.6179 2.679 One T6574 

Bea 28 Ul -0.010 Go ok 1.029 12879 
ez25> 00 Ul 8.389 2-700 O2087 Leoa3 
WIZ Aad 1.940 O2392 12593 

is -3.904 0.898 02056 0.948 

wee 2 Ul le 1.946 0 S6O02 1.3385 
e275 128 Ul =69.839 C2567 0012 027538 
U2 -1.600 PS 450 0.879 32951 

meo 00 Ul -1.679 3.442 0.556 oo 
iz -6.840 36S O20zs 12k? 

B26 12 Ul -0.724 22506 0.834 ones 
226 12B Ul -7.656 SoS OO Ss 1.896 
U2 0.993 4.948 0.854 Ze oo 

Us 5246! 0.828 O7027 0.910 

Ez? 00 Ul -4.194 3,963 0.184 eo 
227 008 Ul 3.048 0.456 0.047 0-675 
WZ 1.424 1.940 0.496 esos 

We =5 7128 2 02076 L.473 

ezei2 Ul ~4.642 5.093 0.192 Zea 
2 0.547 02067 O.214 02285 

eze-00 Ul 122509 5.660 OF ds5 2.379 
U2 8.149 2325 0.034 2525 

U3 2.342 9.411 0.643 3.068 

U4 -4.650 3.155 0.148 Poss 

226 00B Ul 0.918 02625 0.431 0.791 
Zo 12 Ul 4.249 1.663 0.084 LaZzes 
Zoo 12B Ul -4.985 5.404 OSLSo 22825 
0% ole: 325079 0.587 L692 

Zezo 00 Ul eer SEO L763 0.487 Lae 
LEZ =—3,573 Brele9 0.289 22263 

229 12 Uil 9.977 Ti 068 0.100 Bao2s 
e2o9 128 Ul — 2a 4303 0.425 2.074 
U2 2.416 eo a 0.214 eZ56 

U3 e293 0.942 O33:635 Osa al 

yoo 00 Ul Sooo 6 25.670 0.661 52067 
UZ 0.908 4.895 02862 2.2038 

is =3 2.536 14.368 0.566 3.790 

Zoo 12 Ul -3.864 10.684 0.419 3.269 
WZ 0.641 1.494 OnroF Liee22 

Us 3.714 SZ? 0.268 2-242 

yoo 128 Ul 0.878 Goo O54 2855 
Zone 00 - Ul 3.290 15.740 0.600 3.967 
Ze2z500 Ul 4.540 1.136 0.052 12.066 
Zoe eOO0B Ul 3.199 Le2)5 0.106 LOZ 
Peewee Ul -4.410 peg ed 0.084 123530 
iz ZG TOR S26 0.700 S245 


Les 





52 
233 


233 


33 


234 


234 
234 
234 


Z 3. 
235 


ZOO 


23.5 


236 


236 


Z35'°6 
23) 


Z31 
ZS 6 


236 
aso 
240 
240 
241 
ZA 


242 


ZA3 


244 


2B 
00 


12 


26 


00 


0OB 
eZ 
iZB 


00 
00B 


IL 


12B 


00 


eZ 


tZ2B 
00 


O0B 
2 


ZB 
00 
00 
0OB 
00 
00B 


00 


00 


eZ 


U3 
Ua 
OL 
UZ 
U3 
U4 
Wd 
WZ 
EE 
U2 
Ofal 
U2 
U3 
U4 
Ul 
WZ 
Ul 
Walt 
WZ 
U3 
a 
Ul 
WZ 
WG 
al 
U2 
U3 
Wa 
jak 
U2 
U3 
Ul 
Uz 
Ul 
Ul 
UZ 
Ul 
Ud 
U2 
U3 
U4 
Ul 
al 
U2 
U3 
Ul 
U2 
Ul 
Ud 
Ut 
2 
Ul 
U2 
is 
oe 
U2 
Ul 
U2 


Morro rhon fer 


oo 
oD 
. 748 
7518 
2629 
on 0 
-409 
~042 
2320 
nooo 
pos 
226 
-269 
e295 
Wa5 3 
-242 
nO5 
.548 
7o 90 
HSS 
-249 
eZ 3 
22.03 
- 664 
#ea19 
-/96 
Bo eZ 
2056 
20 
mI0Z 
s0S0 
ZUG 
ao 04 
sol 6 
oreo 
ss 
255 
O26 
-248 
.194 
2760 
AY, 
2168 
- 407 
a7 
-462 
aoe 
.341 
-oS8 
Oi 2 
.014 
324 
20.50 
roe 
soo! 
. 768 
oa 
ie Ps 


[= NO 
NO OWA POKRPNWWOWADKHFPHPHWOHND POON NWNHDAWOKRPNBRDITOKPNAHONO OP HAAN SPW BOON S 


W 


— 


Co W 


ke 


fo 


fa 


0 
mos 
-940 
~918 
eos 
O03 
oo 
io 
7520 
so o4 
e210 
.208 
26S 
oe 
ae o2 
By 
- 186 
-064 
-467 
.016 
-461 
a2 i 
wOZZ 
noo 
21:08 
1226 
26) 
TTA 
- 438 
woos 
b255 
-433 
70 
2005 
.920 
2234 
7309 
Ook) 
7911 
fo2) 
.474 
eo 
a3 
eo 
974 
SATAN 
wOns 
eo 
pret 
-496 
a 
oS 
ao 
-986 
ao.o 1 
037 
"og 
860 


LOG 


oooeoercdcooo0oacoaooaoacaoaoaoaaaoaaoaoaorccocorcccdooaoceodocococececececececeo0oncneoooo0o ooo oe 0 9°o oO °° 


lO 
we 9 
BOS 
woe 
a GG 
sou/ 
eS 
~942 
eee 
ey 
2260 
.846 
ailear6 
-449 
082 
a3 
neve) 
i249 
One 
oot 
BON) 
UO 
onlay 
Ss 
~343 
2 
- 348 
2c 
woo8 
nee 
2015 
~402 
- 906 
009 
027 
moUS 
eo 
oles 
eS 
225 
7002 
oa 5 
eee 
-042 
Soo 
790.0 
eo 
50.65 
-489 
oo 
a2ekg 
.144 
Obs 
BO 6)al 
898 
a7 
-440 
oe 


PWHrRNHYONOPHPHENONNNFPWHEWNHRKEPNWORPUOWHN UWOWKHNNNUOPRPNNWNHWNHR BNNRFPNRPFNN WWE DH 


.184 
OT, 
SUieo 
aoe 
a0 60 
so 0 
949 
sS'6 
: 929 
.314 
987 
a2 0 
-401 
-440 
-145 
.274 
.284 
20 LS 
350 
~420 
2209 
2 
cao dl 
oe 
o207 
-492 
2300 
880 
a 
.943 
-495 
2600 
ooo 
-44] 
SoD 
S59 
soo 
Gaal 
02 
.886 
ORL 
bo 20 
.334 
why 
om 
SOG 
oo 
ZOL5 
-864 
600 
oS 
=o ol 
OI 
a2 oS 
-614 
»143 
moe 
oie: 





244 


245 


246 


247 


247 


248 


248 


248 


248 


249 


249 


ZOO 
250 


Zo 1 


Zo. 


Ze 


ZZ 
ZS 


i253 


00 


00 


00 


O0B 


00 


00B 


eZ 


28 


00 


00B 


00B 


2 


12B 


00 


00B 
00 


V Component 


DAY/TIME 


224 


224 


224 
ZO 


22) 
22 


226 


00 


EZ 


12B 
00 


2 
i2B 


00 


vi 
V2 
val 
V2 
V3 
V1 
Vl 
V2 
eS 
V1 
ial 
V2 
val 


cee i, 
70159 
862 
-946 
RO 
.041 
OO) 
22S 
7200 
aod 
-481 
SUZ 
2158 
- 664 
-764 
-092 
woG? 
9:0 
-940 
1250 
noon: 
.974 
Odes 
3290 
.874 
cau 
ous 
aoo4 
soley 
-407 
.614 
5 ly 
ol! 
-485 
760 
oe 
oS 


MEAN 
(m/s) 


mMeRNOWONNOOON 


1 
NO 


oe) 
.174 
eZ 
.776 
-274 
so 
-041 
coe 
-401 
804 
0 G7 
-481 
ah 59 


ra 


So Se Gn cn 


1S, 
ez 


14 


34 
16 
ea 


18 


Bh , 
See 


505 


218 
.006 
-948 
~442 
O62 
-647 
moO 
2a 1 
a Od 
[256 
O16 
OR 
JE(0) 
- 486 
- 806 
47. 
ot 
-008 
Nar 
ee 
.141 
Oy 
=u 2 
~454 
-642 
Boe 
slg) 
vole 
167 
ale 
LEG 
noo 
TOG 
OR 
.670 


889 
913 


829 


445 


419 


430 


148 


822 


VAR 
(maz, Sc) 


PmHNHOrF OF OL HOO OC CO 


.449 
7 U0 
HOD 
-492 
-548 
~945 
sO00 
1208 
ed 
.014 
SO o18 
ers 
oor 


Oy, 


Sooo Ss SSeS S SO: OC CO. UD 6S hL@@ O O 6O hO ULhS OO 2 Oe 2 Oo eo a. eS 


2302 
“960 
~1747 
229.8 
004 
3 0 
.184 
20S 
2029 
.463 
- 984 
050 
1035 
.198 
aes 
ald 
~343 
-482 
aol 
eae 
oo 
- 336 
a3 OF 
no OZ 
«106 
-214 
.043 
047 
eee 
2075 
23 
206 
oT Te 
elas 
2635 
257 
OS 


WGTVAR 


To 2S 2 2 GF co oOo 0 oo © 


“078 
cS 
- 420 
-453 
~474 
-465 
025 
.124 
.419 
> 
~6 20 
seo 
-498 


~994 
O29 
;455 
2003 
.244 
“333 
S56 
3263 
ee 
.042 
2325 
-498 
73 
300 
oS 
eon 
.848 
Lo 
ote 
S153 
“a0 
2024 
.068 
200 
EOD 
S206 
0d 
2852 
039 
woo 
O18 
25 
.084 
.274 
wile 
0 
916 


rPW~JI BrPHO AW DOOR RPWRPWWRrRPRDWRWWR WP WOrRNNFRFF AO 


SDEV 
(m/s) 


ao 0 
Oot 
eed: 
a Or 
ols 2 
i32 
2 7 0 
Sao8 
sot 
007 
43/9 
£85 
es 0 


NNR rRPORPONNODWOCOCO 


eZ 





22.0 
Z20 


22] 


22] 


ZZ] 


228 


Z20 
228 
228 
229 
Z2.0) 
wae 9 
230 


230 


230 
234 
ZZ 
EOYs 
SZ 


Z32 
233 


ZS 


2353 


a3 A 


234 


234 
234 


239 


12 
2B 


00 


00B 


2 


00 


OOB 


LAS} 


00 


re 


25 


00 


eZ 


126 
00 


LZ 


2 


00 


O0B 


Zz 
iZB 


00 


VZ 
Vl 
Va 
V2 
V3 
V1 
Vil 
WZ 
Wo 
Vl 
Nez 
val 
V2 
eS 
V4 
Vil 
V1 
Vil 
V2 
Vil 
V2 
Vl 
Val 
V2 
V3 
Wi 
V2 
V3 
V1 
V2 
Vic 
V1 
V1 
a 
vil 
V1 
V2 
wes 
Wilk 
al 
V2 
V3 
V4 
V1 
V2 
V1 
V2 
Vil 
V2 
vs 
V4 
V1 
V2 
Vi 
V1 
V2 
vs 
Vi 


NORPHONFPRPODANRPRPRPNN BREF NEFNDOWUWHBRP OR ADAN HO 


a03 
mele 2 
=o 05 
oD | 
- 166 
OZ 5 
moo 2 
S235 
- 483 
7 © 
2058 
7 OG 
2208 
Oy) 
aol 
260 
old S 
.041 
2055 
.854 
65 
BG Z 
2007 
7305 
oD 
~421 
-364 
73,010 
20.141 
-640 
DO 
5559 
e290 
e200 
S2AkG 
soc 
SP 
woe 
e912 
D610 
moO 
30 
s0S7 
.264 
sib ope 
O50) 
sd) 
393 
eo 
-862 
Arg 
2050 
27d 
027 
als 
1248 
-403 
088 


NO 


Sooo mn Oo Maa IO oo COW OA he ow Or ON G&A OPN TP RrPoO@awrPrPrrPaorPrPSf& Wo OP OO OM Www eo ht 


NO 


FS 


Fe 


be 


Kk 


sO 5 
794 
-495 
el 
220 3 
nOO7 
2.6 
U0 
027 
-254 
.041 
2905 
yo Z 
oo 
20 
-465 
a2 0S 
008 
226] 
2 1. 
~842 
- 450 
si 54 
200 6 
2526 
ROSS 
OS 
soe 
-250 
20-4 
yo l 
-864 
EG 
od 
661 
-994 
-045 
- 600 
oo 
poze 
soe. 
sols 
oD 
7900 
2575S 
-964 
2525 
pool 
EZ 
-364 
a8 
~943 
565 
OG 
300 
2926 
dg 
Ooo 


ee ee eee Gs Se GG eS) OS Se) SO oO fF oe © Oo oF OOo oO Or Coc CO OrFr OO O° © © 0 oOo @ 200 © ao & 


al 
.006 
-858 
22106 
Oe 
~260 
orn 
RG 
834 
ao? 
1035 
.643 
ot85 
nooo 
ia 
wooo 
-009 
=027 
9 
nS 
Oa 
=o 
Oak9 
foe 
066 
aS 
e038 
~439 
-008 
Ze 
7 0 
SOs 
2092 
Ole 
000 
soo 
eau 
Poo! 
glee 
O22 
SSNS 
s76 
oeleO 
oO 
ooD2 
aS 
1.61: 
SS 
-244 
—oikZ 
3510 
55 
soa 
s505 
= OnE9 
ESIC 
2059 
eZ 


©—OOrFROONRFRRFRF OF 


SDWonmnowWrPnn BNOBOORrRN BENPFPORPORPNNORPRPRPNRPRPRPNRFPRPRPRPRrFORRN 


aod 
O20) 
.870 
co26 
s097 
560 
Zo 
soos 
. OES 
.504 
5203 
-089 
.214 
vouee 
“300 
oe 
1.90 
.004 
ens 
2059 
2260 
2107 
S116 
72510 
ReoZ 
2108 
Zou 
Po” 
woo? 
620) 
20d 
2.05 
-473 
ao OS 
7 9 
sso 
21D 
aol 
~349 
BoSZ 
A eye 
~954 
eo 
-949 
.108 
oe 
EOS 
-206 
«425 
23 
23 09 
- 456 
.O73 
7260 
2309 
A OV 
2206 
aoe 





ZOD 


23.) 


235 


250 


Z306 


Z36 
Ze] 


aS | 
Zo 


ZS 
239 
240 
240 
241 
241 


242 


Ze 3 


244 


244 


245 


246 


247 


247 


248 


248 


248 


248 


249 


00B 


eZ 


P28 


00 


2 


23 
00 


O0B 
EZ 


12B 
00 
00 
00B 
00 
00B 


00 


00 


A 


2B 


00 


00 


00 


O0B 


00 


00B 


dZ 


ZS 


00 


V1 
V2 
V3 
val 
V2 
V3 
wel 
il 
V2 
V3 
V1 
V2 
vil 
V1 
V2 
V1 
V1 
V2 
V3 
v4 
ve 
V1 
V2 
V3 
vl 
V2 
V1 
ee 
v1 
V2 
wal 
V2 
V3 
V1 
V2 
V1 
V2 
V3 
V1 
V2 
V1 
V2 
V3 
V1 
V2 
V3 
val 
V1 
V2 
vl 
V2 
V3 
vl 
V1 
V2 
vl 
V2 
Wak 


{ 
MOnDrOrFPHWUUOEFWPBRRPODRFPODRPRPHRPNFHFNOrF ON 


to 
Pe eae 


MON~ TDOrP~I10OW~INOORFORPRFWNHN OF © 


sey) 
2504 
716 
.014 
od 
-446 
nOZ7 
oO! 
5 (isis: 
2/06 
7 80 
BOGS 
spied 
109 
~445 
oe 
~347 
-461 
7 GZ 
- 924 
W262 
7269 
sE2ro1G 
aoe 
moo 
318 
~477 
20 
-402 
-492 
oo 
e245 
.054 
.004 
-470 
2S 
BOG 
hel 
- 864 
.846 
So 
~145 
eal 
Sone 
5 
-042 
od 
OZ 
S145 
a0 
nOZal 
. 1/64 
biel 
-241 
woo 
OZ 
2096 
20 Go 


me 


NO 


te 


te 


te 


kK 


St 


- i 
S 
Sete ne wt EtG or OP SN OU OOO WN DON TW eo IP ONNOAEP ON N@P hw N Pe ore N 6 oo er icc 


na 


fe 


SR 
ooo 
se Ol 
702 0) 
s00 1 
mod” 
bol 
ZOD 
GZ 
-o87 
eZ 02 
7963 
cZo) 
2108 
19:2 6 
035 
2070 
Bs Ae 
o> 
-247 
oe 
7/786 
- 636 
2219 
-467 
~943 
.026 
2 4 
109 
pono 
-242 
woz 
-491 
oe 
-489 
Sous 
-470 
7220 
.479 
008 
-215 
e250 
-401 
oO 
Lao 
oc 
~/63 
m7 Go 
7.08 
eae) 
~441 
-306 
ooo 
oe 
5 eS) 
ce 
a0 
26Gb 


SSeS ere ea oO. Sm S| OS" O OO @. 2 2 oD © CoCo oOo 0 Oo Cc Oo Oo 0 0 0 O° CC O00 © oO © 0 © 0 0 C0 2 a2 © 2 2 r ® © 


5806 
994 
2429 
O21 
2003 
eoed 
oOo 
7600 
~457 
5169 
.060 
-246 
oe 
27OL 
255 
ooo 
~O50 
3 
786.0 
2916 
2059 
noe, 
wOZe 
21 I6 
~342 
2/76 
oon) 
-841 
08 
948 
wO5I6 
oo 
7200 
-482 
936 
5 ake 
07 
Use 
2662 
—ouS 
poo? 
alibs 
2Ode2 
=0 36 
~947 
20 16 
. 906 
20335 
Bae) it 
Oa 
0'65 
10 18 
OS 
OS 
- 430 
Oo 
007 
ne ou) 


WNHODNNOKFPKFNHNORF OND DCHNHRFONOCWOOCKRrFRKRPHNHNFRBPNFRFNWRONDKHFKNFKRNFRFRRFRRWOPrRPrRPPN POR HN LS 


43 
sao 
304 
e213 
S11) 
O47 
-488 
-456 
wo 7 
ao 
-449 
-459 
eZ 
3603 
sod 
1042 
~945 
208 
2659 
ea 
25009 
605 
od 
910 
SA 
-438 
OS 
. 148 
2D) 
oe 
soo 
607 
-914 
EOS 
2578 
eo 
. 686 
-469 
-934 
ee 
2204 
208 
7550 
eS 
-084 
~ 743 
602 
2256 
oa 
cA 
105 
.818 
2082 
213.6 
7 O63 
~854 
735 
e260 


yz. 
eZ 
eZ 

49 





249 
Z50 
250 


a 1 


Zo i 


Zaz 


OZ, 
Zo 3 


00B 


00B 


PZ 


128 


00 


00B 
00 


Vz 
V1 
V1 
V1 
V2 
Vs 
V1 
vl 
V2 
V1 
V2 
Vil 
Vi 
V2 
VS 
V4 


.019 
.009 
. 768 
soe 
229 
. 070 
314 
-Ouke 
oo. 
188 
036 
40D 
-420 
-J25 
83 
.654 


ss 2 as 
~905 
iO 
010 
.744 
Gy 
207 
2oO9 
odes 
ee 
soa 
lees 
.841 
29S 
2605 
.079 


110 


oo © — ©@ 2 Co eo 2 oo: © 2] 2 & 


One 
.642 
Gs 
oS 
-144 
2) 
mee 
- 900 
69) 
Gi 
-419 
704 
015 
934 
.814 
fot 


rPNHOWOMRPNWOWOaONrRRNN 


. 108 
2OZ.0 
-042 
-418 
oar 
sees 
Os 
0G 7 
7259 
SUEZ 
e208 
OG 
Ot 
me oZ 
Ou 
2039 


81 


ZL 
64 
S56 


V2 
100 
49 
49 
49 
56 


49 
49 
S6 





Table 9: FNOC component wind statistics - 1-day averages 


DAY MEAN VAR WVAR SDEV N 
(m/s) (m*/s*) (m/s) 
224 U -0.234 20.735 1.000 4.954 378 
V 1.193 3.309 0.701 1.819 378 
eon) -1.229 23-510 0.943 4.886 340 
V 0.758 7.354 0.930 2( 12 340 
Zo -2.223 18.383 0.789 4.288 463 
V -0.201 13.628 0.999 3.692 463 
Zot 0.624 ooeZO7 0.991 5.938 499 
V -0.820 9.059 0.933 3.010 499 
228 U 2.234 31.682 0.866 5.629 453 
V -0.013 7.074 1.002 2.660 453 
229 U 0.807 21.933 0.973 4.683 631 
V 2.245 9.824 0.662 3.134 631 
230 U 0.222 20.806 0.999 4.561 695 
V 2.980 7.705 0.465 2.776 695 
231 U 3.772 6.800 0.330 2.608 134 
V 3.080 6.322 0.401 2.0914 134 
one) 1.556 18.062 0.884 4.250 466 
V On 9.150 0.789 3.025 466 
233 U 1.190 16.892 0.924 4.110 636 
V 0.532 11.917 0.978 3.452 636 
234 U -0.568 14.678 0.980 3.831 623 
V 0.455 9.704 0.981 3.115 623 
Zoo U 0.696 20.454 0.978 4.523 662 
V 1.184 10.423 0.883 3.228 662 
236 U 1.452 19.291 0.903 4.392 693 
V 2.038 10.751 Oi 3.279 693 
237 U 0.582 10.881 0.974 3.299 206 
V 0.518 1.926 0.881 1.388 206 
238 U -1.388 18.064 0.905 4.250 945 
V 3.324 6.577 Vie 2.965 945 
239 U -1.528 20.189 0.899 4.493 372 
V 3.513 Tato 0.386 2.184 372 
240 U -2.285 8.209 0.612 2.865 412 
V eT 13.837 0.922 3.720 412 
241 U -2.497 5.200 0.455 2.280 453 
V 0.052 14.181 1.002 3.766 453 
242 U -1.982 6.738 0.633 2.596 260 
V -0.312 9.412 0.994 3.068 260 
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Table 10: FNOC component wind statistics 


1-day /10° latitude averages 





20° — 30° N 
DAY MEAN VAR WVAR SDEV N 
(m/s) (m*/s*) (m/s) 
224 U -4,231 4.351 0.196 2.086 eZ 
V 0.040 0.941 1.004 0.970 172 
D205 -4,965 9.001 0.268 3.000 15g 
V 0.384 T3Ol3 0.986 2.704 esr 
226 U -5.617 5.126 0.140 2.264 234 
V -2.169 G.3a2 0.575 2.516 234 
U2 -3.559 10.310 0.450 a2 11 257 
V -1.502 7.547 Oni 2 2a 257 
228 U -1.197 14.172 0.911 Soo 255 
V 0.258 3.090 0.983 5S 255 
229 U -2.290 11.989 0.698 3.462 250 
V 0.783 3.943 0.868 1.986 255 
230 U -2.846 14.322 0.640 3.785 291 
V 1.306 4.004 0.703 2.001 291 
231 U none 
V none 
Zoe -3.658 3.864 0.224 1.966 134 
V 1.331 3.192 0.646 eon 134 
ao) -2.869 6.129 0.428 2.476 183 
V 1.016 2.349 0.697 1.533 183 
2340) -1.812 8.910 O33 2.985 268 
V 0.421 3.92% 0.960 1.982 268 
235 U -1.686 5.977 0.680 2.445 234 
V 0.302 4.197 0.983 2.049 234 
230 -0.257 4.509 0.989 2123 298 
vy 0.550 2.169 0.880 1.473 298 
Za -0.508 6.168 0.965 2.484 170 
V 0.393 1.603 0.917 1.266 170 
238 U -1.276 8.599 0.844 2.932 234 
V 2.004 5.945 0.598 2.438 234 
239 U -2.051 7.415 0.640 223 206 
V 2.164 4.995 0.517 2230 206 
240 U -3.048 5.489 0.372 2.343 ao 
V -0.115 6.111 1.001 2.472 wo 
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Table 13: Boxed mean differences - 18 Aug 1978 (Day 224) 


Mean Atlas SASS winds - 10 deg lat/10 deg lon boxes and total 
EOredata trom revs 7445. 750) 751, 758 


LONI LON2 LAT1 LAT2 UMEAN UVAR VMEAN VVAR N 
60 70 40 50 2333 Gao oS 0290 ini 
bo 60 40 50 Zo Gao7 Ase A 4.78 2 
40 50 40 50 3205 i ao o de Se oF cle nS 
70 80 3.0 40 25 3.47 Peo2 4.88 Zt 
60 70 30 40 Ze OZ S53 Su 8 S138 89 
50 60 50 40 rr 0 6221 ea 2.6 6:95 94 
40 50 30 40 = Oror, LO Ooo 6.68 81 
70 80 20 30 -4.44 Zn 50 1.00 Sac 49 
60 70 20 30 -4.54 a 2 Zot a5 0. 116 
50 60 20 30 mes) Gro —On a4 4.41 Ze 
40 50 20 30 Oe as This sh), =o 2396 64 
40 80 20 DO Or ll Zio Da OF) 25 6 


Mean FNOC winds - 10 deg lat/10 deg lon boxes and total for 12Z 18 Aug 


LON1 LON2 LAT1 LAT2 UMEAN UVAR VMEAN VVAR N 
70 80 40 50 S242 Om 2 -2.86 23390 Za 
60 70) 40 50 5.42 Ales Gz 6 29.43 it 
DIO 60 40 DO 4.64 Caso 7204 Z oe lide 
40 50 40 S10 S250 G05 ae 14ae05 20 
70 80 30 40 1.64 4.99 0.03 AS Eze 
60 70 30 40 Zag OFZ 0 4.79 4.00 AZ Ak 
D0 60 310 40 1.44 0.66 Zo ZA0.3 v2 
40 50 30 40 05-75 is 164 0.42 2a 
70 80 20 30 -2.45 Lao 7 or Zee Za 
60 20 20 30 -2.75 OR 2.14 3.24 ia 
50 60 20 30 -3.20 O27 0.26 1.49 HZ 
40 50 20 30 —-3.86 10.68 =). 01 oS 2a 
40 80 20 50 0.98 ens 1.47 PSersoe te 7 

Difference of boxed component means: 
Atlas-FNOC: U component/V component 
x8Qx x7TO* x60* x5Q* x4Qx 

x50x 
N/A =3.09/5.34 Jeo fe oo —2,51/0.33 

x40* 
8299/17. 19 0263/0298 -0.04/-1.55 -1.42/-1.05 

x*320* 
-1.99/-0.51 -1.79/0.19 1.83/-0.82 -0.07/-0.64 

x2?0* 
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Table 14: Boxed mean differences - 24 Aug 1978 (Day 230) 


Mean Atlas SASS winds - 10 deg lat/10 deg lon boxes and total 
foOrsdata from revs “830 636, e637, 844 


LON1 LON2 LAT1 LAT2 UMEAN UVAR VMEAN VVAR N 
60 70 40 50 Sic o0 42.43 = 04 57236 6 
50 60 40 50 2 OW 12241 eee tO So 
40 50 40 210 ORG” 36.64 50.3 Pees» 101: 
60 70 30 40 0.54 108.5 Aas 4.35 67 
50 60 30 40 4.86 62.19 2.45 LO 30 65 
40 50 30 40 S70 396 0.45 4o2 i110 
70 80 Z0) 30 Soa Z 0.00 Die 16 0.00 ds 
60 70 20 30 -0.41 Dees 502 Seit 104 
20 60 20 30 =2.42 oa Oat 4.45 ID 
40 50 20 30 =2.7 0 2 G7 eral S220 33 
40 80 ZO 0 ed 6 26.007 161 ats6 ot) OLS 


Mean FNOC winds - 10 deg lat/10 deg lon boxes and total for 122 24 Aug 


LON1 LON2 LATI1 LAT2 UMEAN UVAR VMEAN VVAR N 
70 80 40 50 eG Zs 0.26 10.47 iZi 
60 TO 40 50 OF eZ 23209 -~3.60 36.39 bau 
0 60 40 50 1276 7.66 =4.,07 23.69 2a. 
40 50 40 50 1.90 S2e3 7 6.07 A 7 21 
70 80 30 40 0.14 9.01 =1,18 3.23 al 
60 70 SO 40 0.42 4.47 -1.29 3223 eZ 
50 60 30 40 356 HOe2Z 1.14 22350 iz 
40 50 30 40 6.19 S02 ee 2 0.47 il 
70 80 20 30 -2.99 ee 5 0x39 2,95 Weak 
60 70 20 30 -0.67 2.43 0.64 OZ ez 
50 60 Za) 30 -1.68 4.27 deeeltey Za eek 
40 50 20 30 -3.02 Tao Ou 0.74 eZ 
40 80 20 50 0.64 23 200 O27 too 27a 

Difference of boxed component means: 
Atlas-FNOC: U component/V component 
*BO* x7Q* *60* x5Q* *4Q* 

*50Qx* 
N/A 251 270 oS esos sil, ZS ae! 

*4Qx* 
N/A 0212702 eS Silas Ora /= 1.45 

x30%* 
=0.73/4.717 0.26/2.99 -0 57-036 Neos / 2472 

x2?0x 


za 





Table 15: Boxed mean differences - 5 Sep 1978 (Day 248) 


Mean Atlas SASS winds - 10 deg lat/10 deg lon boxes and total 
foumadedarrom revs 002, 1003," 1006, 5 1009, 1016 


LONI LON2 LATL LAT2 UMEAN UVAR VMEAN VVAR N 
60 70 40 50 6.40 S220 0.34 Lao9 hal 
40 50 40 50 8.66 3926 10.48 26.85 104 
70 80 30 40 Saleen ee Ze —4.9 6.49 he 
60 70 30 40 4.66 4.43 One 4 2 13.64 54 
50 60 30 40 SoZ Gs Gaol 2ss0 34 
40 50 30 40 4.50 4.55 FOS Ga2o. “100 
70 80 20 30 =4 550 4.28 4.64 0.44 8 
60 70 20 30 =2.03 LO; 86 2.49 Seo. LOS 
50 60 Z0 30 =o Ss OS eee eo 
40 50 20 30 =07.41 a9 5 Ss 6.41 53 
40 80 20 50 O59 Ones S200 2929 OU. 


Mean FNOC winds - 10 deg lat/10 deg lon boxes and total for 12Z 5 Sep 


LON1 LON2 LAT1 LATZ2 UMEAN UVAR VMEAN VVAR N 
70 80 40 50 142 19.44 eo ZOroS eZ 
60 70 40 50 4.84 14.65 -2.11 16.94 2a 
50 60 40 a0 4.62 24.93 a2 0 ZeiG9 ea 
40 50 40 50 JasAL9 A283 TO 44.39 iil 
70 80 30 40 =1.06 6770 =0.97 Gavel Za 
60 70 30 40 S220 5.00 = Oe 5eo5 Za 
50 60 30 40 ele. 6:00 S00 Gros 12 
40 BY 30 40 5.60 ere OF 5.51 4.66 pot 
70 80 20 30 eae 0.42 ie 2S 0.90 2 
60 70 20 30 -1.94 Solo 20 On 73 Za 
50 60 20 30 -2.42 TeAG ioe 4.11 leak 
40 50 20 30 oi LOR SS eo: 4.92 eal 
40 80 20 50 oS 29.98 lizv6 5 Zo 8 ee? 1 

Difference of boxed component means: 
Atlas-FNOC: U component/V component 
*B0* x7Q* x*60* x*5Q* x*4Q* 
*5QO* 
N/A 1 5 :Gy2 ea N/A a7 7 2.78 

*40* 

=O O/—5.70 Mest /-0.35 =256672e06 eiloy—-1.44 
*30* 

eee 3. 40 =0.09/=0. 2153 a = 2. 92/1.74 
k20* 
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Table 17: Atlas-FNOC difference field component statistics - 1-day averages 
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DAY MEAN VAR WVAR SDEV N 
(m/s) (m*/s*) (m/s) 
224 U 0.126 5.692 1.000 2.386 378 
V 1.355 6.333 0.777 2.516 378 
225 U -0.917 9.039 0.859 2.245 340 
V 0.686 4.447 0.907 2.109 340 
226 U -0.357 4.217 0.973 2.054 463 
V 0.000 6.518 1.002 2.993 463 
227 U 0.598 5.827 0.944 2.414 499 
V -0.266 5.887 0.990 2.426 499 
228 U 0.695 7.846 0.944 2.801 453 
V -0.913 Sli 7 0.909 2.860 453 
229 U -0.730 9.868 0.950 3.141 631 
IV -1.307 7.141 0.808 2.672 631 
230 U -0.412 11.391 0.987 3.375 695 
V -1.480 5.720 0.724 2.392 695 
231 U -0.233 6.686 0.999 2.586 134 
V -0.766 2.903 0.837 1.704 134 
232 U -0.618 6.341 0.945 ZO ls 466 
V -0.644 10.690 0.965 3.270 466 
233 U 0.065 10.250 1.001 3.201 636 
V -1.410 12.745 0.866 3.570 636 
234 U 0.498 6.919 0.967 2.630 623 
V -1.116 10.804 0.898 3.287 623 
235 U 0.234 5.698 0.992 2.387 662 
V 0.198 9.399 0.997 3.066 662 
236 U 0.275 7.607 0.992 2.158 693 
V -0.283 10.829 0.994 3.291 693 
237 U 1.148 6.725 0.840 2.593 206 
V -2.644 5.662 0.449 2.380 206 
Zao U 0.666 7.995 0.949 2.821 545 
V -0.111 7.885 1.000 2.808 545 
239 U 0.743 6.408 0.923 2.531 372 
V 0.158 3.277 0.995 1.810 372 
240 U -0.564 3.733 0.924 1.932 412 
V 0.486 9.241 0.959 2.289 412 
241 U -0.463 12.716 0.986 3.566 453 
V -0.321 9.003 0.991 3.000 453 
242 U -0.301 19.147 0.999 4.376 260 
V -2.572 8.591 0.566 2.931 260 





243 U 


244 U 


245 U 


246 U 


247 U 


248 U 


249 U 


250 U 


251 U 


252 U 
V 


-1.033 
-3.132 
-0.054 
Sebi 
=1.128 
-2.237 
-1.044 
-0.441 
-1.421 
-0.263 
-0.839 
0.230 

-0.036 
-0.505 
0.484 

-1.211 
1.031 

0.959 

-0.429 
1.318 


2.19% 
7.982 
31.594 
24.144 
4.227 
5.790 
14.779 
15.709 
10.770 
11.738 
15.860 
16.622 
9.672 
7.143 
7.264 
3.925 
26.992 
26.897 
16.367 
15.125 
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0.678 
0.451 
1.003 
0.954 
0.587 
0.537 
0.934 
0.990 
0.843 
0.996 
0.958 
0.998 
1.002 
0.968 
0.972 
0.730 
0.964 
0.969 
0.992 
0.899 


1.482 
2.825 
5.621 
4.914 
2.056 
2.406 
3.844 
3.963 
3.282 
3.426 
3.982 
4.077 
3.110 
2.673 
2.695 
1.981 
9.195 
9.186 
4.046 
3.889 


98 

98 

340 
340 
259 
259 
378 
378 
612 
612 
920 
920 
423 
423 
257 
207 
404 
404 
317 
317 





Table 18: Atlas-FNOC difference field component statistics - 1-day/10° latitude averages 











20° — 30° N 
DAY MEAN VAR WVAR SDEV N 
(m/s) (m*/s*) (m/s) 
224 U -0.834 Zee 0.788 1.589 iz 
V 0.385 2.672 0.953 1.635 ii 
225.U -0.786 2.504 0.806 1.582 157 
V 0.490 2.817 0.927 1.678 157 
226 U -0.597 3.039 0.898 1.743 234 
V 0.042 4.983 1.004 22a 234 
Zod -0.150 5.024 0.999 Da Zot 
V -0.549 4.780 0.944 2.186 257 
228 U O82 8.660 0.993 2.943 255 
V -1.336 7.479 0.810 2130 255 
229 U -1.142 4.845 0.790 2.201 2a5 
V -0.435 5o34 ll 0.969 233i 255 
230 U -0.041 6.681 1.003 2.585 291 
V -0.167 gros 0.996 1.878 291 
Zou none 
Vv none 
232 U -1.392 1.507 0.439 1227 134 
V -0.466 2.299 0.920 1.516 134 
250 U -0.066 3.592 1.004 1.895 183 
V -0.156 4.212 1.000 2.052 183 
234 U 0.262 4.423 0.988 2.103 268 
Vv 0.699 3.986 0.894 1.997 268 
235 U 0.002 6.887 1.004 2.624 234 
V 0.914 4.452 0.845 2.110 234 
236 U 0.852 5.839 0.892 2.416 298 
V -0.832 9.161 0.933 SO 298 
237. 12555 6.933 0.745 2633 170 
V -2.634 5.784 0.456 2.405 170 
238 U 0.017 3.931 1.004 1.983 234 
V 0.298 3.675 0.980 Only 234 
239 U 0.141 3.034 0.998 pig 2 206 
V 0.551 2.529 0.897 1.590 206 
240 U -0.601 3.785 0.915 1.945 Soe 
V 0.816 5.679 0.898 2.383 ont 
241 U -1.202 4.650 0.765 2.156 291 
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242 U 
243 U 
244 U 
245 U 
246 U 
247 U 
248 U 
249 U 
250 U 
251 U 


252 U 


30° — 40° N 


DAY 


224 U 
225 U 
226 U 
227 U 
22510 
229 U 
230 U 


231 U 


0.787 
-2.386 
-0.642 
-1.576 
-1.124 
0.065 
-0.846 
-0.649 
-2.851 
-1.603 
-0.528 
-1.526 
-0.120 
-0.635 
0.138 
-0.018 
-0.662 
0.400 
-1.527 
2.938 
2.239 
0.223 
0.511 





6.242 
3.179 
3.082 
0.371 
2.734 
12.965 
22.692 
3.739 
7.633 
3.584 
8.471 
3.757 
6.578 
8.079 
5.891 
11.852 
8.654 
8.597 
3.852 
23.631 
14.432 
3.268 
5.137 


VAR 
(m*/s*) 





4.891 
8.987 
6.977 
5.559 
4.449 
5.533 
4.304 
4.316 
6.582 
7.590 
7.064 
6.544 
9.006 
6.050 
2.757 
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0.913 
0.360 
0.890 
0.130 
0.694 
1.007 
0.976 
0.905 
0.486 
0.584 
0.971 
0.618 
1.001 
0.955 
1.000 
1.003 
0.955 
0.987 
0.625 
0.735 
0.745 
0.993 
0.964 





0.656 
0.647 
0.965 
0.975 
0.991 
0.953 
0.587 
1.000 
0.643 
0.750 
0.963 
0.564 
0.955 
0.613 
0.911 


2.498 
1.783 
1.755 
0.609 
1.654 
3.601 
4.764 
1.934 
2.763 
1.893 
2.910 
1.938 
2.565 
2.842 
2.427 
3.443 
2.942 
2.932 
1.963 
4.861 
3.799 
1.808 
2.395 


291 
98 

98 

49 

49 

134 
134 
134 
134 
306 
306 
342 
342 
359 
359 
310 
310 
193 
193 
185 
185 
121 
121 


170 
170 
183 
183 
193 
193 
193 
193 
198 
198 
219 
219 
219 
219 
85 





73200 


233 U 


234 U 


235 U 


236 U 


237 U 


238 U 


239 U 


240 U 


241 U 


242 U 


243 U 


244 U 


245 U 


246 U 


247 U 


248 U 


249 U 


250 U 


Zab 


Zon 
V 


-0.289 
-0.468 
-0.128 
0.482 
-1.380 
0.701 
-2.347 
0.410 
-0.092 
0.638 
-0.837 
ns 1 
-3.138 
0.873 
-0.663 
0.915 
-0.557 
-0.422 
-0.783 
-2.281 
-0.651 
-2.123 
-3.011 
-0.490 
-5.141 
1.097 
-0.418 
-2.385 
-2.035 
1.333 
-0.072 
-0.582 
-1.095 
-0.209 
-1.072 
0.810 
0.343 
0.734 
-0.260 
0.790 
-1.442 
-0.277 
1.135 


2.587 
7.248 
7.583 
12.026 
11.087 
9.741 
11.678 
5.157 
6.947 
5.152 
Oa 
5.343 
Soe 
9.268 
12.782 
7.131 
5.386 
3.950 
1.549 
5.902 
1.235 
6.622 
7.764 
3.466 
5.162 
66.644 
37.485 
1.197 
3.735 
56.002 
46.855 
Sento 
30.133 
7.590 
6.211 
6.214 
3.039 
3.231 
2.985 
18.939 
19.662 
26.826 
22.043 


135 


0.980 
0.974 
1.002 
0.984 
0.855 
0.955 
0.681 
0.975 
1.002 
0.929 
0.932 
0.631 
0.275 
O9Zi 
0.970 
0.904 
0.956 
0.963 
0.723 
0.517 
0.750 
0.599 
0.463 
0.953 
0.164 
0.990 
1.004 
0.174 
0.477 
0.982 
1.014 
0.996 
0.968 
0.997 
0.846 
0.910 
0.969 
0.869 
0.993 
0.978 
0.913 
1.004 
0.951 


1.608 
2.692 
2.754 
3.468 
3.330 
321 
3.417 
2.399 
2.636 
2.270 
3.028 
2.311 
1.929 
3.044 
3.575 
2.670 
2.321 
1.884 
1.245 
2.346 
Peel 
2.573 
2.186 
1.862 
2.212 
8.164 
6.123 
1.094 
1.933 
7.483 
6.845 
5.726 
5.489 
2.159 
2.492 
2.493 
1.743 
1.797 
1.728 
4.352 
4.434 
5.179 
4.695 


85 
268 
268 
3953 
393 
270 
270 
294 
294 
338 
338 
157 
157 
279 
279 
85 
85 
85 
85 
98 
98 
98 
98 
49 
49 
121 
121 
eZ 
(2 
12 
t2 
149 
149 
355 
355 
149 
149 
64 
64 
98 
98 
147 
147 





40° — 50° N 


DAY 


224 U 


225 U 


226 U 


227 U 


228 U 


229 U 


230 U 


Zale) 


232 U 


233 U 


234 U 


235 U 


236 U 


Zan U 


238 U 


239 U 


240 U 


241 U 


242 U 


243 U 





VAR 


(m*/s*) 


4.131 
3.175 


5.718 
1.533 
7.798 
16.747 
7.798 
16.747 
18.249 
8.480 
16.196 
5.273 
13.213 
2.414 
17.260 
32.106 
15.074 
23.696 
5.881 
9.209 
5.375 
15.493 
12.262 
14.743 


13.359 
7.306 
11.638 
2.060 


32.619 
8.524 
32.619 
8.524 
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0.446 
0.484 


0.588 
0.090 
1.021 
0.985 
1.021 
0.985 
0.847 
0.816 
0.873 
Oral 
1.013 
0.492 
0.953 
0.983 
0.892 
0.874 
0.973 
0.600 
0.856 
1.002 
1.000 
1.002 


0.954 
1.003 
0733 
1.008 


0.936 
0.467 
0.936 
0.467 


SDEV 
(m/s) 


2.032 
Loz 


2.391 
1.238 
292 
4.092 
2.792 
4.092 
4.272 
2.912 
4.024 
2.296 
3.635 
1.554 
4.155 
9.666 
3.883 
4.868 
2.425 
3.035 
2.318 
3.936 
3.502 
3.840 


3.655 
2.103 
3.411 
1.435 


5.711 
2.920 
5.711 
2.920 


36 
36 


36 
36 
49 
49 
49 
49 
157 
157 
234 
234 
49 
49 
100 
100 
172 
172 
121 
121 
234 
234 
242 
242 


162 
162 
81 
81 


113 
113 
113 
113 





244 U 


245 U 


246 U 


247 U 


248 U 


249 U 


250 U 


251 U 


Zou, U 


none 
-1.880 
-2.517 
-3.716 
-0.855 
none 
none 
-2.156 
0.356 
-1.404 
2.487 
none 
none 
none 
none 
-1.689 
0.946 
-2.499 
3.865 


6.471 
5.160 
2.304 
0.900 


ea 
2.613 
39.855 
38.426 


29.996 
45.078 
12.473 
9.818 
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0.652 
0.451 
0.143 
0.558 


0.330 
0.961 
0.957 
0.864 


0.908 
0.989 
0.676 
0.400 


2.044 
ea lee 
1.518 
0.949 


1.509 
1.616 
6.313 
6.199 


9.099 
6.714 
3.532 
3.133 


85 
85 
49 
49 


121 
121 
242 
242 


121 
121 
49 
49 





Table 19: Atlas-FNOC difference field vector statistics - area averages 


REV TIME MEAN VARSPD SDSPD VARVEC SDVEG N 
p28 1 224 00 S174 4.997 ZS Zone 31560 49 
Z 1760'> OR G52 Ors 923 Zoe eronno 36 
664 1 224 12 Zoe 5 {76 3 2.157 T3629 Ze Oe 36 
Z et) eo eee lees 3.144 ed eS 64 
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Table 23: FNOC wind stress curl - individual area averages 
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Radar geometry for scatter from a sea surface, 


(from Stewart, 1985) 





Figure 3: 
NORTH POLE 
Ww 
<, 
“Ad 
a 
—~ 
= 
O 
z 
3 EAST 
see DOIRBLUER GELL 
ys MEASURE MENT 
; LOCATION 
ae oe PROJECTION OF ANTENNA LOOK 
DIRECTION ONTO EARTH SURFACE 
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Figure 5: Examples of noise-free solutions for beam patterns 90° 
apart, vertical polarization, and a wind speed of 10 m/s 


for wind directions of X = 0°, 40°, and 80° 
(from Boggs, 1982) 
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Figure 10: Distribution of ungridded Atlas wind components - averaged 
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Figure 12: Distribution of mapped Atlas wind components - averaged 
over individual boxes. a. U velocity, b. U stnd dev., 
c. V velocity, d. V stnd dev. 
N = 145 pts 
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Figure 15: Mapped vs. ungridded Atlas wind component standard 
deviations - from l-day averages 
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Figure 16: Mapped Atlas vs FNOC wind speed andmcdtmectilome—- 1-day perizod 
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Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - 1-day period 
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Figure 19: Mapped Atlas vs FNOC wind speed and direction ~ 1-day period 
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Figure 20: Mapped Atlas vs FNOC wind speed and direction - 1-day period 
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Mapped Atlas vs FNOC wind speed and direction - 1-day period 
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Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Figure 23: Mapped Atlas vs FNOC wind speed and direction - 1-day period 
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Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - 1-day period 
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Figure 28: Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - 1-day period 


180 





ATLAS - FNOC WIND COMPARISON 
DAs 9 e372) PTs 


FNOC WIND SPEED (M/5) 





0 5 10 15 20 
ATLAS WINO SPEED (M/35) 


240 


180 


120 


60 


FNOC WIND DIRECTION 





0 } 60 L120 180 240 300 360 
ATLAS WIND DIRECTION 


Figure 31: Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - 1-day period 
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Figure 33: Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - i-day period 
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Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - |-day period 
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Mapped Atlas vs FNOC wind speed and direction - |l-day period 
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Mapped Atlas vs FNOC wind speed and direction - I-day period 
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Mapped Atlas vs FNOC wind speed and direction - 1-day period 
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Mapped Atlas vs FNOC wind speed and direction - l-day period 
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Mapped Atlas vs FNOC wind speed and direction - 1-day period 
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Figure 45: Ungridded Atlas winds - Rev 1008 
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Figure 46: FNOC winds - Day 248 122 
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Figure 47: Atlas vs FNOC vector speed - l-day averages 
N = 29pts 
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Figure 48: Atlas vs FNOC mean speed and vector standard deviations 
I-day averages N = 29 pts 
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Figure 50: Atlas vs FNOC mean speed and vector standard deviations 
3-day averages N = 10 pts 
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Figure 52: Atlas vs FNOC component standard deviations 
l-day averages N = 29 pts 
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Figure 53: Atlas vs FNOC component velocity - l-day/10° lat averages 
20°N - 30°N N = 28 pts 
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Figure 54: Atlas vs FNOC component standard deviations 
l-day/10°lat averages 20°N - 30°N 
N = 28 pts 
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Figure 56: Atlas vs FNOC component standard deviations 
l-day/10° lat averages 30°N - 40°N 
N = 29 pts 
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Figure 58: Atlas vs FNOC component standard deviations 
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Figure 59: Distribution of all individual Atlas-FNOC difference field 
vectors. a. AU velocity, b. AV velocity, c. vector speed, 
d, vector direction 
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Figure 60: Atlas wind stress curl - Rev 664 boxes 2 and 3 
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Figure 61: Atlas wind stress curl - Rev 715 1 
(x 10? N/m? ) 
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Figure 62: Atlas and FNOC wind stress curl - Rev 1059 2 
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Figure 63: Atlas and FNOC wind stress curl - Rev 830 1 
Gaon Ny mo 


ZS 





He Cee aes a1 0-6 


— 
: 
a 


\ 
47 N ce ‘ 
* \ \ w~ 
Osa on MED 
oF Nias 
45°N o 
Sh 
= 4 
2 45'N a 
= o> 
c ; 
— 44° O 
ae — ° cae 
: ‘O 
4356N wn 
2.5 Ra 
a: 
42 N aS 
41° N 4 4 3 4 a a 4 4 4 
49 Ww 48 W 47 W 45 WW 45 ww 44° WN 43 W 42 W 4] WW 
LONGI TUDE 
1008 2 FNOC xX 10-65 
OS 
q..0 
43'.N “OSD § 
-1.0 | 
47".N i 
» 
46° N -1.9 > 


LAT I TUDE 

&> ‘ 

e S 
9°24 
(\ 

© 

an) 

O 





1.0 
a = 
“— on 
‘ Oo 
N+ : 
42 en 
41 N a 4 A 4 a 1 4 9 4 
49 WN 48 W 47 W 45 45 WwW 44 WwW 43 OW 42 W 41 «# 


LONGITUQe 


Figure 64: Atlas and FNOC wind stress curl - Rev 1008 2 
( x 10°* N/m? ) 
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Figure 65: Atlas and FNOC wind stress curl - Rev 1008 1 
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Figure 66: Bragg scatter from short waves tilted by longer 
waves in the composite surface approximation, 
The vertical scale is exaggerated. 
(from Stewart, 1985) 
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